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X XV. BAND SECHSTES (SCHLUSS-)HEFT 1957 


Experimentelle Durchfiihrung einer neuen hydrodynamischen 
Analogie fiir das Torsionsproblem * 


Von G. Gro®mann 


1. Einleitung. Eine Reihe wichtiger Probleme der Elastizitatstheorie laBt sich mathematisch 
nur mit hohem Aufwand lésen. Hierzu gehért auch die Bestimmung des Spannungszustandes eines 
auf Torsion beanspruchten prismatischen Kérpers. Die Torsionsspannung lat sich nur in wenigen 
Fallen berechnen, so da schon friih — vor etwa 50 Jahren — Gleichnisse entwickelt wurden, um 
auf experimentellem Wege zum Ergebnis zu kommen. Die bekanntesten Analogien sind das 
Membran-Gleichnis (Seifenhaut) von Prandtl und das Thomsonsche Strémungsgleichnis. Vor einigen 
Jahren veréffentlichte Pestel’ eine neue hydrodynamische Analogie, deren experimentelle Durch- 
fiihrung im folgenden beschrieben wird. 

In seiner Arbeit betrachtet Pestel die Bewegung einer zihen Fliissigkeit zwischen einer ebenen 
Scheibe und einer hierzu parallelen Ebene. Bewegt sich die Scheibe unter der Einwirkung einer 
Kraft P langs einer Normalen auf die Ebene zu, so wird die Fliissigkeit aus dem Spalt h zwischen 
der Scheibe und der Ebene herausgedrickt. 

Die Differentialgleichungen dieser Bewegung liefern folgende Beziehung zwischen dem Druck p 
unter der Flache und der im Spalt h zwischen Scheibe und Ebene durchflieBenden Fliissigkeits- 
menge q: 

ope 12 Op 127 


Ox a h3 Wx? By = AB dy a (1) 


Diese Gleichungen lassen sich unter Benutzung der Kontinuitatsgleichung wie folgt zusammen- 
fassen: 
e Ge 12 
P+ pba h. (2) 
(eee ake h® 


Dabei ist 7 die dynamische Zahigkeit der Fliissigkeit. 
Vergleicht man (1) und (2) mit den Differentialgleichungen fiir die Torsion prismatischer Stabe 
nach de St. Venant 


Op Op 

Bat ey? By tex? (3) 
a. eee 
Ox2 | By RES (4) 


so kann man auf Grund der Gleichheit des Aufbaues die folgenden Zuordnungen erkennen: 


Torsion Strémung } 
69>. 
—26%9 ——+> aa ht (a), 
i$ ——> (»). 
GGin¥ + iy dF =M, > fpdF=P (), | (5) 
F F 
6 
Ty m Ee “4 Vx (d) 7 
6 
oe, a ne — = dy (e) ¢ J 


* Auszug Dissertation Hannover 1957; Referent: Prof. Dr.-Ing. E. Pestel, Korreferent: Prof. Dr.-Ing. 
E. Martyrer. Die Arbeit entstand im Festigkeitslaboratorium der TH Hannover. 
1 E. Pestel, Ing.-Arch. 23 (1955), S. 172. 


bo 
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2, Aufstellung der Gleichungen fiir die experimentelle Bestimmung des Torsionswiderstandes. 
Bei einem prismatischen Kérper der Querschnittsflache F’, der sich unter Wirkung der Kraft ¢ 
in einer Fliissigkeit vom spezifischen Gewicht y und der dynamischen Zahigkeit 7 mit der Ge- 
schwindigkeit h° bewegt, greifen unter Vernachlassigung von Reibungskraften an den Seiten- 
flachen des Zylinders die in Abb. 1 eingezeichneten Krafte an. Nach dem d’Alembertschen Prinzip 
JaBt sich daraus folgende Differentialgleichung aufstellen: 


—_mh'+ P+A—G=0. (6) 


Darin bedeuten m die Gesamtmasse, die an der Bewegung teilnimmt, G die Kraft, die die Bewegung 
erzeugt, A = y F (H —h) den Auftrieb und P die durch die strémende Flissigkeit hervorgerufene 
Kraft. Wenn man bedenkt, daB 2 G ) = 2 M,/Jp gilt, worin Jp der Torsionswiderstand des unter- 
suchten Prismas ist, lat sich fiir P aus (5a) und (5c) der folgende Ausdruck ableiten: 


Dae (7) 


Mit A und P in (6) eingesetzt lautet die Gleichung 


EES sf the (6 (8) 


h? 


mh” 


Fiir diese Differentialgleichung existiert keine geschlossene Lésung. Sie kann deshalb nur numerisch 
aufgelést werden. 

Eine Voraussetzung fiir (5) ist das Vorhandensein schleichender Strémung, d. h. die auftreten- 
den Fliissigkeitsbeschleunigungen miissen so klein bleiben, daB die Massenkrafte vernachlassigt 
werden kénnen. Nach Pestel ist diese Bedingung er- 
fullt, wenn 


Ret een 


v 


ist. Die numerische Auflésung von (8) bestatigte 
dieses; denn nach einer von der Zahigkeit 7 des Oles 
und dem Torsionswiderstand Jp der Flache F' abhan- 
gigen Anlaufzeit wurde das Beschleunigungsglied so 
klein, da es vernachlassigt werden konnte. Glei- 
chung (8) geht dann itiber in 


Abb. 1. so, GS Gey ie 
— Sie i (9) 


Fir den Fall, daf8 man den Auftrieb vernachlassigen kann, nimmt (9) die einfache Form an 


G 
ies Mere eth a 
h 3 Ip ae : (9a) 


die mit Gleichung (7) tibereinstimmt, wenn G durch P ersetzt wird. In dieser Form erméglicht sie 
eine schnelle Beurteilung der Versuche, da sie bei Auftragung der Geschwindigkeit h’ iiber der 
Abszisse h im doppelt-logarithmischen Mafstabe Geraden gleicher Neigung tg x = 0,333.. ergibt. 
Der Abstand dieser Geraden ist durch den Parameter G/3 Jp 7 gegeben. 

Kinfacher als die Messung der Geschwindigkeit als Funktion des Weges ist die Weg-Zeit-Messung. 
Es empfiehlt sich daher, Gleichung (9a) iiber die Zeit zu integrieren. Daraus erhalt man 


1 1 2G 


hi hg ~ FJp-y 1%) - or 


Diese Form 1aBt sich mit den Koordinaten t und h im doppelt-logarithmischen MaBstab auch als 
Gerade auftragen, deren Giiltigkeitsbereich durch die Forderung nach schleichender Strémung 
gegeben ist. 

Fiir den Fall, daB die Wirkung des Auftriebes nicht vernachlassigt werden darf, kann man 
Gleichung (9) in die folgende Form umschreiben: 


_yFht PB 


ale 3Jpy °° 39pn (9b) 
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Darin ist P) = G —y H F die Differenz zwischen der Kraft G und dem Auftrieb A bei verschwin- 
dender Spalthéhe h. Die Integration von Gleichung (9b) ergibt 


1 1 2yF/(1 =| yF\) h(hyF + P) 
Ui i) Py (i hy : ( Py 2 hy (hy y F + Po)I eo 
In dieser Gleichung kann der dritte Term in der eckigen Klammer vernachlassigt werden, da er 
ungefahr eine Zehnerpotenz kleiner als der vorhergehende Term ist. Damit lautet Gleichung (10b) 


anal Le ta Fe si 
Salers as 


In dieser Form wurde sie der experimentellen Bestimmung des Torsionswiderstandes zugrunde 
gelegt. 


3 JD 


3. Aufstellung der Gleichungen fiir die experimentelle Bestimmung der Torsionsspannung. Die 
Gleichungen (1), (5d) und (5e) sind die Ausgangsgleichungen zur Bestimmung der Torsionsspannung. 
Es muf daher entweder das Druckfeld oder die Fliissigkeitsmenge q nach GréBe und Richtung 
an jeder Stelle unter der Platte ermittelt werden. Eine Messung des Druckes p wahrend des Be- 
wegungsvorganges ware experimentell nur mit sehr groBem Aufwand und wenig Aussicht auf 
Genauigkeit méglich gewesen. 

Nun ist zu bedenken, daB die Spannungsverteilung lediglich eine Funktion der Querschnitts- 
form ist, oder anders ausgedriickt, daB der Verlauf der Spannungstrajektorien lediglich eine Funk- 
tion der Querschnittsberandung ist. Sie sind eindeutig 


bestimmt, wenn die zu ihnen senkrecht verlaufenden t 

Linien bekannt sind, deren Richtungen im (x, y)-Ko- 0 

ordinatensystem durch — Tae Ty gegeben ist. Die Glei- ~— 
_ Z aS} 


chungen (5d) und (5e) zeigen, daB diese Richtungen 
gleich q,/q, sind und mit den Stromlinienrichtungen 
des Strémungsgleichnisses zusammenfallen. 

Dadurch bietet sich als weiterer Lésungsweg die Abb. 2. 
Bestimmung des Stromliniennetzes an, welches — 
ebenso wie das Feld der Spannungstrajektorien — nur eine Funktion der Kérperform ist. Die 
hiermit entstandene Méglichkeit, sich bei der Spannungsbestimmung mit Hilfe der Analogie auf 
die Untersuchung der Kinematik des Strémungsvorganges zu beschranken, wurde durch den 
Umstand begiinstigt, daB der absolute Betrag der Geschwindigkeit den Stromlinienverlauf nicht 
beeinfluBt. Werden diese bei sehr kleinen Geschwindigkeiten der Platte mit Aluminiumpulver 
sichtbar gemacht, erhalt man, wenn die Geschwindigkeit im Laufe der Bewegung bis auf Null 
herabgesetzt wird, ein stehendes Bild des Liniennetzes, das sich auswerten laBt. 

Es ist zweckmaBig, an Stelle des («, y)-Koordinatensystems ein krummliniges (s, t) zu wahlen, 
bei dem die Richtung s mit der Richtung der Stromlinien zusammenfallt. Dann ist wegen der 
Definition der Stromlinien q, = 0 und die ,,Torsionsspannung* 7 der Analogie 


af 


a 67 
eras Ror iar (11) 


Hierbei ist g, die Fliissigkeitsmenge, die in der Zeiteinheit in Stromlinienrichtung durch die Lange I, 


in Spannungsrichtung gemessen, flieBt. Die Kontinuitatsgleichung liefert entsprechend Abb. 2 


PO ad Ey 
Lee o.di 


Setzt man hier fiir h’ den Wert aus (7) ein, so erhalt man 


CHE JT? 


poe E 12 
1s ~ di Spas ue) 
Die Gleichungen (11) und (12) ergeben dann 
= P dk ‘ 
= 2 13 
2“ Jp dl (13) 
Wir definieren nun als bezogene Modellspannung t* 
s + JD dF 
St 14 
ee pp 2: (14) 
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Bei einem Kreisquerschnitt wiirde damit die eee Modellspannung am Aufenrand 


Vergleicht man diesen Wert mit den wirklich beim Kreisquerschnitt auftretenden bezogenen 
Spannungen am Au®enrand 


s Jb 

ta = Te Mr 

so sieht man sofort, daB die bezogenen Werte vollig ibereinstimmen. Dieses gilt nicht nur fiir den 

Kreisquerschnitt und nicht nur fiir den apereand sondern fiir jeden Ort i beliebigen Quer- 
schnittsflachen. Also ist 


ea Rs 


ay tea /D. (15) 


Hier zeigt sich die Einfachheit der Analogie. Zur Bestimmung der Torsionsspannung eines beliebig 
peionnten Prismas wird das Stromlinienbild aufgenommen, nach Gleichung (14) die bezogene Mo- 
dellspannung ermittelt and aus dieser nee Glei- 
chung (15) die wirkliche Spannung errechnet. 

Da es fiir die Lage der Stromlinien gleichgiiltig 
ist, ob die Platte sich nach oben oder nach unten 
bewegt, wurde hier im Gegensatz zur Vorgehensart 
bei der Bestimmung des Torsionswiderstandes der 
Versuchskérper gehoben; denn so lieBen sich die 
Stromlinien besser sichtbar machen. 


4. Beschreibung des Versuchsaufbaues und der 
Durchfiihrung der Versuche. a) Versuchsein- 
richtung. Die Versuchsapparatur besteht aus 
zwei nebeneinanderliegenden Wannen aus Plexi- 
glas. Beide haben einen Einsatz aus Messingblech, 
der in der Hohe verstellbar ist, so daB diese 
Platten unabhangig voneinander und unabhangig 
von den Wannen waagerecht ausgerichtet werden 
kénnen (Abb. 3). Bei der Bestimmung des Tor- 
sionswiderstandes iiber die Messung der Fallzeit 

ADU. oo Vedeatlsapmermiar von oben ceschan: dient die Messingplatte gleichzeitig als Grund- 
kontakt. Der Versuchskérper ist titber eine Schnur 
mit einer Waagschale verbunden, die-es erméglicht, beliebige Werte von P, [vgl. (10b)] her- 
zustellen. Eine zweite Schnur fiihrt zu dem Registriergerat, das auf einen Streifen gleichzeitig 
den Weg und die Zeitmarkierungen aufzeichnet. Eine weitere Fiihrung ist nicht erforderlich, da 
der Sinkvorgang stabil ist. Die Ubersetzung des Weges betrug 4:1, d. h. die Héhe wurde vier- 
fach tiberhéht vom Schreibgeraét aufgezeichnet, um den Bereich kleiner Héhen besser auswerten 
zu kénnen. 

Vom MeBschrieb laBt sich die Kurve auf doppelt-logarithmisches Papier iibertragen. Dieses 
ist besonders bei kleinen Werten von h vorteilhaft, da hierdurch die Lange der Kurve reduziert 
wird, Man kann dann sofort iibersehen, ob die Messung mit der Theorie iibereinstimmt (Abb. 5). 

Ein an einem verstellbaren Stativ befestigter Fotoapparat erméglicht es, die einzelnen Versuche, 
besonders aber das Stromlinienbild zur Spannungsbestimmung, bildlich festzuhalten. Versuchs- 
flissigkeiten sind Cylanto 220 fiir die Widerstandsbestimmung und ESSTIC 50 fiir die Herstellung 
des Sitomlunenbilden Beide Ole wurden feed hoeecene von der ESSO A.G. Hamburg zur 
Verfiigung gestellt. Zwei U-Rohrmanometer gestatten es, wahrend der Versuche mit Helles 
sahalarsa en Innendruck zu messen. 

Die gesamte Versuchsapparatur ist auf einem Tisch aufgebaut und mit diesem fest verbunden. 

b) Versuchskérper. Die Versuche werden mit 12 verschiedenen Querschnittsformen durch- 
gefithrt. Diese hatten die in Tabelle 1 angefiihrten Abmessungen. 

Alle Versuchskérper waren aus shrrga elas Pheper Plexiglas heroes tellt bei den Formen 2, 4, 6, 
8, 9, 10 und 11 war auf die Flache ein 12 mm bzw. 15 mm fee und 12 mm hoher igen von 
Plexiglas aufgeleimt worden. Jeder Kérper hat an seiner Bodenflache drei 1 mm lange Stahl- 
spitzen. Sie Ahonen zur Begrenzung der Sinkzeit, die sonst nach Gleichung (10c) gegen mend 
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a 


gehen wiirde. Gleichzeitig bilden sie Grundkontakte, die das Aufsetzen auf die Messingplatte auf 
dem MefBschrieb markieren. 


Tabelle 1. (D = Durchmesser, s = Kantenlange, 6 = Wandstdrke) 


— 
in) 


ee 
FPOoUOmMDAAMNEWwWhNre 


Kreis 1 
Kreis 2 
Quadrat 1 
Quadrat 2 
Rechteck 1 
Rechteck 2: 
Dreieck 1 
Dreieck 2 
Dreieck 3 
Dreieck 4 
Dreieck 5 
Spiralbohrer 


D 
Da 
Ss 


Al TS 
oe 
5 
5 


150 mm 
150 mm 
133 mm 


i 
Se 


Dg = 150 mm 


omaem 


i tl ll 


—~o 


Oo = 15 mm (Hohlquerschnitt) 


12 mm (Hohlquerschnitt) 
04 mm. 
04 mm 6= 12 mm (Hohlquerschnitt) 


12 mm (Hohlquerschnitt) 
28 mm (Hohlquerschnitt) 

5 mm (Hohlquerschnitt) 
12 mm (Hohlquerschnitt) 
. Abb. 9 und 10) 


oe 
ish, 


Die Probekérper sind im allgemeinen an drei, die beiden Quadrate an vier Punkten aufgehanet. 
Die einzelnen Aufhangungen sind in ihrer Lange veranderlich, so daB vor den Versuchen die Grund- 
flache des Versuchskérpers parallel zum Boden ausgerichtet werden kann. 

c) Bestimmung des Torsionswiderstandes (Versuchsergebnisse). In Tabelle 2 sind die 
Ergebnisse der Versuche zur Bestimmung des Torsionswiderstandes wiedergegeben. Zum Vergleich 
wurden die errechneten Werte fiir J) mit angegeben. Fiir die Falle, fiir die keine exakten Gleichun- 
gen zur Berechnung des Torsionswiderstandes vorlagen, wurde dieser naherungsweise nach den 
Bredtschen Formeln ermittelt. Als dritter Wert ist in Tabelle 2 der giinstigste Bereich fiir hy an- 
gegeben. Die Kenntnis dieses Bereiches ist wichtig. Die Versuche zeigten namlich, daB der Anlauf- 
vorgang bei den auftretenden Versuchszeiten einen so kleinen EinfluB auf die Gesamtzeit hat, 
da8 er praktisch unberiicksichtigt bleiben kann, wenn hy sich in den angegebenen Grenzen bewegt. 

Der Versuchskérper hat dann bei h ~ 2 mm immer die Grenzgerade erreicht, die durch Glei- 
chung (10a) gegeben ist, wenn dort hy = oo gesetzt wird. Das bedeutet, daf} die Zeit vom Versuchs- 
beginn, d. h. dem*Loslassen des Kérpers, bis zum Versuchsende im Rahmen der MeBgenauigkeit 


lorsionswiderstand Ip 


°C O/temperartur toe 
seks falize/t te 
10° 


Belastung des Versuchskorgers be 
3 5 678914 


i 
7 


te 4h 


2F-7-7 [mm] 


Band | ani 
ae 


Z 


ia Z 
a yaa 
oS 


Cs! 


ay 


ES 


ay 2 Ty 
04 


£ 


Beispiel: Dreieckiger Hohlzylinder 


& 


S,=dbamm, O=Z8mm 
toe=18,4°C 


_ 


We 
SS 


A,=3009 | 
bf = 160 Sek 
(Nur giiltig fur Cylanto 220) 


‘Abb. 4. Bestimmung des Torsionswiderstandes aus MeBwerten. Nur fiir Cylanto 220 giiltig. 
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unabhangig von der Anfangshéhe hy ist. Auf Grund dieser Versuchsergebnisse 1aBt sich die Glei- 


mm 4 /allhone h 
7 ——_—1— 1 2s 
c + ghee (Zur Auswertung nicht herangezogen) ee eae 
8 : 0-8 fe 4 He at 
7 Wi, Mg=6 
LMig=5 i} 
6 — Ze Fo 0 
; TS 
y 2 Grenzgerade nach Gleichung (10d) 
is === + 
; OS tn 88.580h Tp-589-10° rom =gored 
i: *) Quest: goma* 
j 1 Reichenbticher: no107d"__| 
Sc. 
Sy 
| es, /alvell ty 
1 2 Pe PAGE. 2 3 & Sek 


Abb. 5. Bestimmung des Torsionswiderstandes bei einem Spiralbohrer. 
Cylanto 220, tg¢ = 19,7° C, Py = 250 g. 


Tabelle 2 


Jp [mm‘] Torsionswiderstand 


Giinstiger Bereich 


gemessen errechnet fir hy [mm] 
1 Kreis | 4,78 -10° 4,97 -10? 4—11 
2 Kreis 2 2,95 -107 2293) LO? 4— 9 
3 Quadrat | 4,30 -10? 4,38 +10? 4—12 
4 Quadrat 2 2d OmmealO 2,2) 107 3— 6 
D Rechteck 1 1,69 -10? 173) 107 4—10 
6 Rechteck 2 13528 10% 1,28 ~ 107 4— 6 
7 Dreieck 1 Sexo, a UM 3,60 -107 4—11 
8 Dreieck 2 2,033 = 107 SOO 5S— 7 
9 Dreieck 3 3,10 -107 3,06 -107 3— 7 
10 Dreieck 4 0,632 - 107 0,611 - 107 3— 7 
11 Dreieck 5 0,572 - 10? 0,538 - 107 3— 6 
12 Spiralbohrer | 0,589-107 | 3— 8 


naked ie SS 


chung (10c) noch etwas weiter ver- 
einfachen: 


edna ees 
es eP (i P, fg) ee 


Fiir die Auswertung dieser Gleichung 
wurde das in Abb.4 dargestellte No- 
mogramm entwickelt. 

Bei mehreren Hohlquerschnitten 
wurde wahrend des Sinkens der In- 
nendruck gemessen; er blieb prak- 
tisch wahrend der gesamten Ver- 
suchszeit konstant. 

Bei den Versuchskérpern 1, 2, 3,5 
und 7, bei denen der Torsionswider- 
stand noch exakt berechnet werden 
konnte, lieferte die Messung um 1 bis 
3 v.H. zu kleine Werte. Bei den 
iibrigen Kérpern waren die Versuchs- 
werte durchweg gréBer als die errech- 
neten. Dieses ist auch verstandlich, 
denn die Wandstarke war hier bereits 
so groB, daB sie nicht mehr von einer 
,,Rohre“ konstanten Schubflusses 
ausgefillt werden konnte. Die Bredt- 
sche Formel basiert aber gerade auf 
dieser Annahme. Der Spiralbchrer 
konnte mit Versuchswerten von 


Quest!§ (Jp = 0,513 10? mm*) und 


! H. Quest, Ing.-Arch. 4 (1933), S. 510, 


 errechnete kurve 


Abb. 6, Spannungsverlauf am Rand des gleichseitigen Dreiecks (oben) und Stromlinienverlauf beim gleichseitigen Dreieck. 
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Reichenbdcher (Jp = 0,544 - 107 mm) verglichen werden. Diese Werte wurden frither im Festig- 
keitslaboratorium der Technischen Hochschule Hannover durch experimentelle Auswertung des 
Seifenhautgleichnisses gewonnen, Abb. 5 zeigt die Weg-Zeit-Kurven des Spiralbohrers. 

d) Bestimmung der Torsionsspannung (Versuchsergebnisse). Beim Kreisquerschnitt 
wurde erwartungsgem4B ein vél- 
lig radialer Verlaut der Strom- omtz- 5 
linien erzielt. Dieses Ergebnis 
ist zwar eine Bestatigung der 
Theorie, sagt aber nur wenig iiber 
die Genauigkeit des Verfahrens 
aus. Dafir lieferte die Unter- 
suchung des Dreiecks 1 bereits 
ein scharferes Kriterium. Aus 
den Stromlinien (Abb. 6) wurde 
eine Kurve der bezogenen Span- 
nung T* ermittelt und mit be- 
rechneten Werten  verglichen. 
Die Versuchswerte stimmten, 
wie Abb. 6 zeigt, mit den be- 
rechneten Werten gut iiberein. 
Dieses war auch beim Quadrat 1 
der Fall, wo die gréBte Spannung 
nach der Analogie T* = 9,05 cm 
gegeniiber dem _ rechnerischen 
Wert t* =8,95cm nur um 
1,l v.H. abwich. Zum Vergleich 
wurden bei diesem Modell (Abb.7) 
die Spannungstrajektorien ein- 
gezeichnet und fiir die auBere 
Stromréhre die Spannung nach 


der Bredtschen Gleichung NX ominien 
a \S 
a XS ee 
18 Jp = 2 Mp = 
p 0(s) 
ermittelt und in Abb. 7 einge- ale 
tragen. 
1 H. Reichenbdcher, Ing.-Arch, 7 


(1936), Seams Abb. 7. Spannungs-, Stromlinien- und Potentiallinienverlauf beim Quadrat. 


\ 
mittlere Spannung in x 
a | der Guberen,Schubrohre” ‘ 


Abb. 8. Stromlinienbild. Quadrat 2. Abb. 9. Stromlinienbild. Spiralbohrer. 
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Die Auswertung der tibrigen Querschnittsformen geschah in gleicher Weise. In Abb. 8 sind 
2. B. die Stromlinien des Quadrates 2 (Hohlquerschnitt) wiedergegeben. Abb. 9 zeigt die Strom- 
linien des Spiralbohrers, Abb. 10 die Auswertung dieses Stromlinienbildes. Die gréBte Spannung 
betrug hier Tmax = 0,666 r G @.  (Reichenbicher ermittelte Tmax = 0,636rG 0, Quest Tmax = 


0,687 GB.) 


be— 75cm al 


Abb. 10. Spannungs- und Stromlinienverlauf an einem Spiralbohrer. 


5. Zusammenfassung. In der vorliegenden Arbeit wurde die experimentelle Durchfiihrung 
des hydrodynamischen Torsionsgleichnisses von Pestel beschrieben. Es beruht auf der Gleichheit 
der Differentialgleichung fiir die schleichende Strémung und fiir die de St. Venantsche Torsion. 
Der Vorteil dieser Analogie besteht einmal darin, daB man Torsionswiderstand und Torsions- 
spannungen unabhangig fiir sich, jeweils in der dafiir giinstigsten Art und Weise bestimmen kann. 
Zum anderen ist der Aufwand an Gerdten und an Rechnungsarbeit sehr gering. Fiir die Aus- 
wertung der experimentellen Ergebnisse bei der Ermittlung des Torsionswiderstandes wurde ein 
Diagramm entwickelt. 

Sowohl bei der Torsionswiderstandsbestimmung wie auch bei der Spannungsermittlung wurde 
eine fiir technische Zwecke ausreichende Genauigkeit erzielt. Vergleichsmessungen an Flachen 
mit bekannten KenngréBen (Jp, tT), die beim Prandtlschen Gleichnis z. B. in der Reichenbécher- 
Apparatur noch nétig sind, eriibrigen sich in diesem Fall. 

Der Torsionswiderstand ]48t sich aus einer Zeitmessung und einer Temperaturmessung er- 
mitteln, die Spannung geht auf die Bestimmung der Flache zwischen zwei Stromlinien zuriick. 


Einfache Umrechnungsgleichungen stellen die Verbindung zwischen Modell und Wirklichkeit her. 
(Eingegangen am 23. Dezember 1956.) 


Anschrift des Verfassers: Dr.-Ing. G. GroBmann, Vlotho (Weser), Bonneberger Weg 37. 
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Uber eine besondere Klasse selbsterregter Schwingungen 
Von K. Klotter und E. Kreyszig 


1. Einleitung und Problemstellung. Das Klassische Musterbeispiel eines Systems, das selbst- 
erregte Schwingungen ausfiihrt, ist der sogenannte van der Polsche Schwinger, dessen Ver- 
halten durch die ,,van der Polsche Differentialgleichung“ 


q—0q(l—o% q?) + 2%q =0 (1.1) 
beschrieben wird. 

Aus der Differentialgleichung (1.1) lassen sich die grundlegenden EKigenschaften dieses Schwin- 
gers in einfacher Weise ablesen: Fiir kleine Werte q ist das zweite Glied (,,. Dampfungsglied“‘) 
negativ; dem System wird Energie zugefiihrt. Fiir groBe Werte von q ist das genannte Glied 
positiv; dem System wird Energie entzogen. Man kann demnach erwarten, daB die sich einstellen- 
den Schwingungen im Laufe der Zeit einem stationdren Grenzzustand, einer Grenzschwingung 
(einem ,,Grenzzykel*‘) zustreben. 

Diese groben Uberlegungen lassen sich erstens verfeinern und zweitens auf allgemeinere Diffe- 
rentialgleichungen, von der Form 


q+ qs(q.q) +#f(q) =9, (1.2) 
ausdehnen. In der Tat bildet das Problem der Existenz und der Eindeutigkeit von Grenzschwin- 
gungen bei Systemen, die einer Gleichung von der Form (1.2) geniigen, den Gegenstand zahlreicher 
und ausgedehnter Untersuchungen, von denen ein groBer Teil dem Sonderfall (1.1) gewidmet ist. 
Derartige Betrachtungen haben das Ziel, Bedingungen fiir die beiden Funktionen f(q) und g(q, q) 
anzugeben, die die Existenz — und gegebenenfalls auch die EKindeutigkeit — der Grenzschwin- 
gungen sichern. Einen leicht zuganglichen kurzen Bericht hieriiber findet man in den Anhangen III 
und IV des Buches von J. J. Stoker!, wo vor allem Ergebnisse von N. Levinson und O. K. Smith? 
behandelt werden; eine weiterreichende Darstellung enthalt das Buch von E. A. Coddington und 
N. Levinson®. 

Es ist klar, daf Existenz- und Eindeutigkeitsbeweise allein keine oder nur ganz unzureichende 
Hilfsmittel an die Hand geben, um das Aufschaukeln oder Abklingen der Schwingung zum Grenz- 
zustand hin im einzelnen zu verfolgen. Andererseits sind exakte Lésungen fiir Gleichungen vom 
Typ (1.2), die eine solche Aufklarung geben kénnten, bisher nicht bekannt geworden. Selbst der 
klassische, einfache Sonderfall (1.1) der Gleichung (1.2) ist einer exakten Liésung nicht zuganglich. 
Will man das Verhalten der Lésungen der Gleichung (1.1) oder anderer Sonderfalle der Gleichung (1.2) 
genauer untersuchen, so ist man auf die Verwendung von Naherungsverfahren angewiesen. 
Zu diesen Verfahren zahlen z. B. die verschiedenen Anwendungsformen der Stérungsrechnung, 
weiterhin Verfahren, die auf dem Gedanken der Variation der Parameter aufbauen und geeignete 
Vernachlassigungen in sich schliefSien (wie etwa die Methode von van der Pol oder die von Kryloff- 
Bogoliuboff) sowie die (meist graphisch verwertbaren) Hilfsmittel, die die Phasenebene anbietet. 
Eine kurze Ubersicht iiber derartige Naherungsverfahren findet man z. B. in einem Bericht von 
K. Klotter*. 

Angesichts der soeben geschilderten Lage ist es von Interesse, besondere Typen von (1.2) zu 
untersuchen, deren Lésungsverhalten analytisch gut erfaBbar ist. Die Differentialgleichung 


i—(cen) £0 —v8 QF + 2h =0 mit p>, (1.3) 


auf die K. Klotter® hingewiesen hat, ist von dieser Art. Sie beschreibt einen Schwingungsvorgang, 
der den durch (1.1) beschriebenen Schwingungen sehr dhnlich ist: Fiir kleine Werte q wird dem 
betreffenden System Energie zugefiihrt, wahrend fiir groBe Werte q Energie entzogen wird. Wenn 


1 J. J. Stoker, Non-linear vibrations. New York 1950. 

2 N. Levinson und O. K. Smith, Duke Math. Journ. 9 (1942) S. 382. 

@ E. A. Coddington und N. Levinson, Theory of ordinary differential equations. New York 1955. _ 

4 K. Klotter, Neuere Methoden und Ergebnisse auf dem Gebiet nichtlinearer Schwingungen, VDI-Berichte, 
Bd. 4 (1955) S. 35. 

5 K. Klotter, Journ. Appl. Mech. 22 (1955) S. 493. 
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auch der Mechanismus des Energieaustausches in seinen Einzelheiten bei beiden Schwingern jeweils 
verschieden ist, so fallt das Bild des Schwingungsvorganges, im groBen betrachtet, in beiden Fallen 
recht abnlich aus. Nach einem Vorschlage von K. Klottert wollen wir (1.3) deshalb als eine ,,.mo- 
difizierte van der Polsche Differentialgleichung“ bezeichnen. 

Es ist bemerkenswert, daB die Modifikation nicht nur in der Abanderung des Mechanismus 
fiir den Energieaustausch besteht, sondern auch darin, daf nun beliebige (und nicht nur lineare) 
Riickstellkrafte zugelassen sind, da an die Stelle von q im letzten Gliede von (1.1) in (1.3) eine be- 
liebige Funktion f(q) tritt. 

Da die Differentialgleichung (1.3) einer relativ einfachen analytischen Behandlung zuganglich 
ist, beruht auf dem Umstand, da® diese Gleichung durch Einfiihrung der neuen unabhangigen 


: ; Pe lixd bares : : 
Veranderlichen V = q? und unter Benutazung der Identitat ¢ = — in eine (abschnittsweise) 


2 dq 
lineare Differentialgleichung iibergefiihrt werden kann; man erhalt namlich 
dV : 
oY — (eg j) B (1 a8 g®) V +28 f(q) =0. (1.4) 


Die vorliegende Arbeit ist der Behandlung dieser Differentialgleichung gewidmet. 

Das Integral von (1.4), das zugleich ein erstes Integral von (1.3) ist, und das zweite Integral 
von (1.3) lassen sich, von trivialen (und praktisch uninteressanten) Sonderfallen abgesehen, nicht 
mehr durch eine endliche Anzahl bekannter Funktionen darstellen. Es ist zu erwarten, daB die 
nachstehend entwickelten Methoden zur Diskussion des ersten und des zweiten Integrals auch in 
anderen Fallen niitzlich sein kénnen, die einer ,,geschlossenen“ Integration nicht mehr zuginglich 
sind. Abschatzungen, wie sie im folgenden durchgefiihrt werden, sind oft den tiblichen Naherungs- 
verfahren zur Lésung gewoéhnlicher Differentialgleichungen vorzuziehen, insbesondere dann, wenn 
es sich darum handelt, einen allgemeinen Uberblick iiber den Verlauf eines Schwingungsvorganges 
zu erhalten. 

K. Klotter hat in der schon erwahnten Arbeit? die Differentialgleichung 


j—Geng) oP + 2h) =0 (1.5) 


auf ahnliche Weise untersucht. Wahrend, mathematisch gesehen, (1.5) ein Sonderfall von (1.3) 
ist (der dem Parameterwert « = 0 entspricht), sind die durch (1.5) bzw. (1.3) (mit « ~ 0) be- 
schriebenen Schwingungssysteme physikalisch durchaus verschieden: Die Differentialgleichung (1.5) 
beschreibt ein System, dem dauernd Energie entzogen wird; es handelt sich bei ihren Lésungen 
also einfach um eine Klasse freier gedampfter Schwingungen, die asymptotisch nach null gehen. 
Wir wenden uns nun der Betrachtung der Differentialgleichung (1.3) baw. (1.4) zu. Eine 
Zusammenfassung der wesentlichsten Ergebnisse findet sich im SchluBabschnitt 13. 


2. Kurze Diskussion der Differentialgleichung (1.4). Indem wir 2 x? durch y* ersetzen, schrei- 
ben wir die Differentialgleichung (1.4) in der folgenden Weise: 


V'—B (sgn g) (L—o @)V+yfq@)=0 mit B>0, (2.1) 


wobei also der Strich die Ableitung nach der unabhangigen Veranderlichen q bezeichnet. 
Lassen wir die Bewegung etwa bei einem Ausschlag q = Q, > 0 mit der Anfangsgeschwindig- 
keit q = 0 beginnen, so wird g zunichst negativ, und wir erhalten aus (2.1) 


Vi +p (l—atg’) V+ y*f(q) =0. (2.2) 


Diese Differentialgleichung beschreibt die Bewegung, solange die Geschwindigkeit negativ ist. 
Aber auch bei positiver Geschwindigkeit ergeben sich keine prinzipiellen Anderungen: Durch Dar- 
stellung der Funktion f(q) als Summe einer ungeraden Funktion f,(q) und einer geraden Funk- 


tion f2()> , 

fi =f, rite ? 
und Einfithrung der neuen Variablen ¢ = — q erhalten wir aus (2.1) fiir den Fall positiver Ge- 
schwindigkeit die Differentialgleichung 


Vit B(L—e?G)V + (f,—f) =0, (2.2’) 


1S. FuBnote 5S. 389. 
2S. FuBnote 5 S. 389. 
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in der alle Funktionen von 7 abhangen und der Strich die Ableitung nach dieser Variablen be- 
zeichnet. Gleichung (2.2’) unterscheidet sich nicht wesentlich von (2.2). Ist f(q) insbesondere eine 
ungerade Funktion, wie z. B. in dem einfachsten denkbaren Fall der linearen Rickstellkraft, so 
sind die beiden Differentialgleichungen iiberhaupt miteinander identisch. 

Unter den Funktionen f(q) sind im besonderen diejenigen physikalisch interessant, die ,,echten“ 
Riickstellkraften entsprechen, d.h. die die Bedingungen 


(a) f(qq)290 fir q>0, 
(b) fq) <0 fir qg<0, (2.3) 
(c) f(q) += 0 

erfiillen. 


Die meisten Aussagen, die im weiteren Verlauf der Untersuchung gemacht werden, beziehen 
sich auf echte Riickstellkrafte, also auf Funktionen f(q), die den Beziehungen (2.3) geniigen, ohne 
daf dies jedesmal angegeben wird. Wo ausnahmsweise von den Beschrankungen (2.3a) und (2.3b) 
abgesehen wird (wie etwa im zweiten Teil von Abschnitt 7), ist dies ausdriicklich bemerkt. 

Kin Wort mag auch am Platze sein iiber die Dimensionen der in (2.1) auftretenden GréBen. 
Deutet man die Differentialgleichung etwa mechanisch!, so stellt q im allgemeinen eine Linge oder 
einen Winkel (im Bogenmaf) dar. Demgema8 erhalt man eine der beiden Spalten in der folgenden 
Liste von Dimensionen: 


Entweder Oder 
dimg = L dimq = 1 
abn, | es ILE I RS ita) A ie 
Bienes == she Ginn og == Il 
dim pr L-% | gins) 5s) 
Mita al ny | tren ave tae hess 
duay 0 dim jf 1 


Um die Gleichung von den Dimensionen zu befreien, wahlt man zwei BezugsgréBen, q, und V,, 
und definiert 


a 7 ee ee ies pee 2 1 (aa 
q=-—; Sy a=a4q, B=BQ, P=yY V,? ae 


Nach Multiplikation mit q)/V, geht (2.1) in eine Gleichung iiber, die anstelle von q, V, a, 6, y, f 
die entsprechenden iiberstrichenen (also dimensionslosen) Gréfen enthalt. La8t man dann die 
Striche der Einfachheit halber wieder weg, so kann man jetzt (2.1) als Gleichung in dimensions- 
losen GréBen ansehen. Von dieser Freiheit werden wir spater (Abschn. 6, 9, 10), wo von Beispielen 
die Rede sein wird, Gebrauch machen, ohne auf die Frage der Dimensionen noch einmal zuriick- 


zukommen. 


3. Die ,,Amplitudenbeziehung“. Den Maximalausschlag des Abschnitts einer Schwingung vom 
Ausschlag g = 0 bis zum nachsten Durchgang durch qg = 0 bezeichnen wir (in etwas ,,liberaler‘ 
Sprechweise) kurz als die (positive oder negative) Amplitude dieses Abschnitts (richtiger ware: 


Schwingungsweite). 
Die Differentialgleichung (2.2) hat die allgemeine Lésung 
y? 
pa EG k) , 3 
Vg) =—F (1) +) (3.1) 
wobei 
q 
T(q) = f s(w) f(w) dw (3.2) 
0 
ist und iy 
s(w) =exp(Bw—cw®) mit c= a (3:3) 


gilt; k ist eine Integrationskonstante. Das Integral I(q) hangt selbstredend von der jeweiligen 
Form von f(q) ab. Ist f(g) ein Polynom, so existiert die Lésung V fiir alle endlichen Werte von q 


und ist beschrankt. 


1 Beschreibt (2.1) einen Vorgang in einem elektrischen Stromkreis, so daf q eine Stromstarke oder eine 
Spannung bezeichnet, so kann man eine Befreiung von Dimensionen auf eine analoge Weise erreichen. 
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Fiir die von einer Ruhelage q = Q, > 0 beginnende Bewegung erhalten wir aus (3.1) die par- 


tikulare Lésung 


52 
V(q) =—~ [10,) — 1(@)] : 3.14 
@) = % (a) — 1) (3.1) 
Diese Lésung gilt, solange q < 0 ist. Wird ein Wert q = Q, erreicht, fiir den V(q) erstmals ver- 


schwindet, fiir den also die Beziehung ' 
I(Q,) eae I(Q;) (3.4) 


erfiillt ist, so ist von diesem Punkte ab bis zur nachsten Nullstelle von V(q) die Lésung der Diffe- 
rentialgleichung (2.2') zu verwenden, die man aus (3.1) erhalt, indem man q durch q und f durch 
fu —feg exsetat. Aus (2.3) folgt, daB Q, < 0 sein mu. 

Mit Hilfe der ,,Amplitudenbeziehung (3.4), die zur Differentialgleichung (2.2) gehort, und 
einer entsprechenden, die zur Differentialgleichung (2.2') gehért, falls (2.2') nicht mit (2.2) identisch 
(also f(g) nicht ungerade) ist, kann man die Folge der Umkehrpunkte ermitteln. 

Um die weiteren Darlegungen einfach zu halten, werden wir explizit nur den Fall ungerader f(q) 
behandeln. Dann sind die Differentialgleichungen (2.2’) und (2.2) und damit die aus ihnen folgenden 
Amplitudenbeziehungen jeweils identisch. Im Falle nicht ungerader Funktionen f(q) miBte man 
die beiden Folgen Q,, Q3,Q;,... und Q5, Qy, Q,,... getrennt betrachten. Das Vorgehen ist in 
diesem Fall jedoch (nach dem, was hier gesagt werden wird) so selbstverstandlich, da wir von 
einer Besprechung im einzelnen absehen kénnen. 

Da eine Schwingung durch die Folge der Amplitudenwerte Q; weitgehend beschrieben wird, 
so ist die Amplitudenbeziehung (3.4) wichtig. Wir widmen ihr daher besondere Aufmerksam- 
keit. 


4, Untere Schranken fiir die ,,Grenzamplitude“ Q*. Wie einleitend begriindet wurde, existieren 
fiir die durch die Differentialgleichung (1.4) bzw. (2.1) beherrschte Schwingung zwei Gebiete, die 
zu unterscheiden sind, namlich ein Gebiet der Energiezufuhr, das bei kleinen Werten des absoluten 
Betrages des Ausschlages q liegt, und ein Gebiet des Energieentzuges, das groBen Werten des 
absoluten Betrages von q entspricht. Die Nullstelle ¢ = |a|-1 des Koeffizienten des Dampfungs- 
gliedes von (2.1) bildet die Grenze zwischen den beiden genannten Gebieten. 

Im Abschnitt 1 haben wir schon plausibel gemacht, daB die Folge der Betrage der Ampli- 
tuden Q; einem Grenzwert Q* zustrebt; wir wollen diesen Grenzwert als die ,,Grenzamplitude“ 
bezeichnen. Wenn (fiir nicht ungerade f(q)) zwei getrennte Folgen zu erwarten sind, so strebt 
jede fiir sich einem Grenzwert zu. Wir unterdriicken jedoch — wie oben gesagt — die Bespre- 
chung dieses Falles, weil er keine weiteren Schwierigkeiten bietet. 

Der Beweis fiir die Existenz eines endlichen Grenzwertes Q* folgt im nachsten Abschnitt. Hier 
wollen wir zunachst eine untere Schranke fiir Q* angeben. 

Zuniachst iiberlegt man sich leicht, daB stets 


OP Saigin (4.1) 
sein muB, falls Q* existiert. Ware namlich Q* < |x|, so wiirde die betreffende Schwingung aus- 
schlieBlich in dem Energiezufubhrgebiet verlaufen, und ein derartiger Wert Q* kénnte nicht auf- 
recht erhalten werden. 

Wir wollen nun eine scharfere untere Schranke fiir Q* herleiten. Hierzu betrachten wir die 
Amplitudenbeziehung (3.4). Die Nullstellen der Funktion 6 w — c w® liegen bei 0 und + wy», wobei 


Wo — -Ve= oss 


| | 


ist. Die durch (3.3) definierte Funktion s(w) erfillt demnach fiir — oo < w < —w, und fir 
0 < w < w die Beziehung 
(Wiad 


In den restlichen beiden Intervallen der w-Achse gilt 
O=<'s(w)iss is 


Auf Grund der Voraussetzung, daB f(q) eine ungerade Funktion ist, erhalten wir damit wegen 
(2.3) die Ungleichung 


I(— 4) < 1(q), (0 <a a (4.2) 


la 
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Beginnt also die Schwingung bei einem positiven Werte Q,, der nicht gréBer als /3/|x| ist, so ist 
die nachste Amplitude, Q,, dem Betrage nach stets gréBer als Q,. 

Demnach ergibt sich: Bei einer Schwingung, die der Differentialgleichung (2.1) 
gehorcht, gilt fiir die Grenzamplitude stets die Abschatzung 


ON se, ibs (4.3) 


Es ist bemerkenswert — und aus der Differentialgleichung (2.1) verstandlich — daB die durch 
(4.3) gegebene Schranke nicht von £ abhangt. Weiterhin ersieht man aus (4.3), daB sich die Grenz- 
amplitude Q* mit abnehmenden Werten von |x| nach gréBeren Werten verschieben muB; auch 
dies ist ohne weiteres einleuchtend. Strebt « gegen null, so muB Q* gegen unendlich streben, d. h. 
es existiert dann keine endliche Grenzamplitude mehr. 


5. Die Beschranktheit der Amplituden Q; und die Existenz der Grenzamplitude Q*. Die Am- 
plituden Q,, Q3,... einer Schwingung, die der Differentialgleichung (2.1) mit 
einem Polynom f(q) entspricht, sind (fiir jede noch so groB gewahlte Anfangsaus- 
lenkung Q,) beschrankt. 

Diese Aussage ist, vom Standpunkt der Mechanik gesehen, fast selbstverstandlich; sie laBt sich 
in einfacher Weise auch formal bestatigen. Ist namlich f(q) ein Polynom, so gilt wegen (2.3) 


lim f(q) s(q)=— oo, lim f(q) sq) =9. 


Gir 5.00 Laat 


Und zwar wachst die Funktion f(q) s(q) fiir groBe negative Werte von q dem Betrage nach sehr 
stark an, wahrend sie fiir wachsende positive gq sehr rasch gegen null strebt. Hieraus folgt unmittel- 
bar die Richtigkeit der obigen Aussage fiir endliche Werte von Q,. Zu zeigen ist noch, daB I(q), 
wenn f(q) ein Polynom ist, auch fiir beliebig wachsendes positives q beschrankt bleibt. Von einem 
gewissen, von f und c abhangigen Werte q =a > 0 an ist B g—cq? < —q. Da weiterhin [(q) 
fiir endliche Werte |q| beschrankt ist, so existiert eine (von f,c¢ und f(q) abhangige) Konstante 
M (< oo) derart, daB 


I(a)< M + fe fw) dw 
gilt. Indem wir f(q) in der Form : 
f(a) = & am q” ee) 


schreiben, erhalten wir demnach 
q n n 
lim I(q)< M+ lim fe’ Ya, w"dw=M+4+ 3 a,0(m+1)<0, (5.2) 
q-> 00 qIr7>~X 90 m=0 m=0 
was zu zeigen war. 

Auf Grund von (2.3) ist I(q) fiir 0 < q < oo und auch fir 0 = q > — oo eine monotone Funk- 
tion von q. Hieraus ergibt sich: Im Falle einer ungeraden Rickstellkraft ist die Folge 
der Betrage der Amplituden |Q,|, |Q,|,..., monoton. é 

Damit erhalten wir zusammenfassend das Ergebnis: Im Falle einer ungeraden ee 
stellkraft existiert fiir jeden Wert x40 eine (endliche) Grenzamplitude Q j 

Ferner gilt demnach fiir zwei aus verschiedenen positiven Ruhelagen beginnende Schwingungen: 


mit 0, > 0, ist auch |Q,| > [2 fiir alle Werte j = 2,3,... 


6. Weitere Eigenschaften der Amplituden Q;. Es ist klar, daB die einer gewissen festen Anfangs- 
amplitude Q, entsprechende nachste Amplitude von der GréBe der in der Differentialgleichung (2:4) 
enthaltenen Konstanten abhangen muB, also von a, / und den in f(q) enthaltenen Koeffizienten. 
DaB der Wert von y belanglos ist, sieht man aus (3.2) und (3.4). Wir wollen die genannte Abhangig- 
keit nun etwas naher untersuchen. 

Fiihren wir in (2.2) die neue Variable 


«=f q (6.1) 


ein, so erhalten wir im Falle einer der Funktion 


ioc} (6.2) 
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entsprechenden Riickstellkraft die Differentialgleichung 
—— (1 = win?) Vi? Bo Ae ee 0 te (6.3) 


Fiir eine von der Ruhelage x = 8 Q, = X, > 0 aus beginnende Schwingung ergibt sich die Sonder- 
lésung 


Via) yn oe aI) ie) (6.4) 


wobei 


De) eae ee ree ears (6.5) 
0 


ist. Die (3.4) entsprechende Amplitudenbeziehung lautet nun: 
J(X1) = J(X;) - (6.6) 


Da das Integral J(x) auBer dem festen Parameter m nur den Parameter u = «/f enthalt, also bei 
Vorgabe eines gewissen Wertes von u wohlbestimmt ist, so hangt das Amplitudenverhaltnis 


IX! _ BGI _ [ul 
[X. ~ FO! ~ [Ql 


nur von u ab. Entsprechend schlieBt man fiir die weiteren Amplituden. — Es ist 


z(u, w) = exp (w — u? w?/3) 


fiir jeden festen Wert w > 0 (bzw. w < 0) eine mit wachsendem u monoton fallende (bzw. wach- 
sende) Funktion. Demnach nimmt J(x) bei monoton wachsendem u fiir jeden (festen) positiven 
Wert von x monoton ab und fiir jeden (festen) negativen Wert von x monoton zu. 

Ks gilt also: Im Falle einer Riickstellkraft, die einer Funktion von der Form (6.2) 
entspricht, ist das Verhaltnis |Q,|/|Q,,1| der Betrage zweier aufeinanderfolgender 
Amplituden nur von dem Parameter , 

Len (6.7) 
abhangig und wachst monoton mit wachsendem u. 

Die GréBe u bezeichnen wir als den wesentlichen Parameter eines durch die Differential- 

gleichung (2.1) beschriebenen Schwingers. 


7. Obere Schranken fiir die Grenzamplitude Q*. Wie gezeigt wurde, existiert fiir « ~ 0 und 
ungerades f(q) eine (endliche) Grenzamplitude, und die Funktion I(q), vgl. (3.2), ist monoton. 

Im folgenden benétigen wir Schranken fiir I(q). Wir geben solche fiir den Fall an, daB f(q) ein 
Polynom von der Form (5.1) mit nichtnegativen Koeffizienten ist. Die im Integranden von I(q) 
auftretende Funktion exp (§ w) wird durch ihren kleinsten bzw. gréBten Wert im Integrations- 
intervall ersetzt. Es ist 

q- 


— cw ».m La) m m+ 1 
[e w dw—=-¢ ( the P Feet) : (7B) 
6 
wobei 
PGn,*) a) (cea 2 oem (7.2) 
0 
die unvollstandige Gammafunktion bezeichnet. So ergibt sich unmittelbar die Abschatzung 
A(q) < I(q) <e? Aq) fir q>0 (7.3) 
mit 
leas - = m+1 
A(q) Sat HS ma * an P 3 seat). 


Wir geben nun obere Schranken fiir die Grenzamplitude Q* an. Hierzu betrachten wir die 
Amplitudenbeziehung (3.4), und zwar der Einfachheit halber fiir die spezielle ungerade Funktion 


f@=@"+t! (m=0,1,...). (7.4) 


\\ 
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(Die nun folgende Uberlegung 1a8t sich auch fiir andere Riickstellkrafte durchfiihren; jedoch 
ergeben sich dann kompliziertere F ormeln.) Wir bestimmen einerseits eine obere Schranke K(q) 
von I(q), q > 0, und andererseits eine untere Schranke L(q) von I(—q). Dann gilt 


I(q) — I(— q) < K(q)—L@) fir q>0. (7.5) 


Wegen (4.2) ist die linke Seite von (7.5) fiir hinreichend kleine Werte q(> 0) stets positiv. Demnach 
ist auch die rechte Seite von (7.5) fiir solche Werte q von null verschieden. Bestimmt man die 
Schranken K(q) und L(q) so vorsichtig, daB einerseits K(q) fiir groBe q langsam und monoton 


anwachst und andererseits L(q) fiir groBe q ein rasches monotones Wachstum besitzt — beides 
ist auf Grund des Wachstumsverhaltens von I(q) und I(—q) méglich —, so hat die Gleichung 
K(q) — Lq) = 0 (7.6) 


fig poe Werte von q(> 0) genau eine Nullstelle, und diese Nullstelle ist cine obere Schranke 
DEO”. 
Wegen (7.3) und (7.4) ist 


uf 
I(q) < 3 cote NS Pp G (m+ 1),¢ ¢) ed = fir g>0 
und 
1 2 
I(—q) > 2m +98 ate acts —eg') e-84  fir—q> 0. 
Also wird 


1 
Iq) — 1(—q) <3 2 +9 (82 P(3 (m +1), 48) —e-Pt P(Z(m +1), —eg')}. 
Die Gleichung (7.6) hat demnach im vorliegenden Falle die Gestalt 
2a P(z(+0),¢9) = P(z(m +1), eg"), (7.6') 


und die Nullstelle g = q) dieser Gleichung ist auf Grund der vorstehenden Uberlegung eine obere 
Schranke fiir die Grenzamplitude Q* des betrachteten Schwingers. 
Gleichung (7.6’) laBt sich durch Anwendung der Funktionalbeziehung der unvollstandigen 
Gammafunktion, 
P(m + 1, x) =m P(m, x) — e—* x™, 
einfach handhaben. Wegen 
Pla) Leer” 


braucht man insgesamt nur die Werte der P-Funktion fiir die Parameter m = 1/3 und m = 2/3 zu 
kennen, um (7.6’) auszuwerten. Diese Auswertung geschieht am einfachsten durch graphische Er- 
mittlung des Schnittpunktes der Kurven, die den Funktionen auf den beiden Seiten von (7.6’) 
entsprechen. Wenn nétig, verbessert man diesen Wert mit Hilfe eines der iiblichen Naherungs- 
verfahren. 

Fiir den soeben behandelten Fall (und ahnliche Falle) lassen sich auch solche oberen Schranken 
von Q* finden, bei denen man nur elementare Funktionen zu betrachten braucht. Dies soll nun 
gezeigt werden. 

Die Funktion f(q) habe die Form (7.4). Wir bestimmen wiederum eine obere Schranke von I(q) 
fiir q > 0, und eine untere von I(— q) und verfahren dann weiter wie oben. Wegen (2.3) und (5.2) 
gilt im vorliegenden Fall 

I(q) < lim I(q)< 0 fir q>0, 
q —> 0c 
und wir kénnen sogar eine von q unabhangige obere Schranke fiir I(q) angeben. Daf eine derartige 
Abschatzung relativ roh ist, wird sich auf das Endergebnis nicht sonderlich auswirken, da I (—q) 
fiir wachsende Werte von q(> 0) sehr rasch anwachst. 
Wegen (7.4) ist 
q 
I(q) = f s(w) we™tldw. 
0 
In geeigneten Intervallen der w-Achse kann die Funktion s(w) bei den vorzunehmenden Abschat- 


zungen durch die Funktion 
Twp) 626” 
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ersetzt werden. Fiir w = 0 und w = -—- W, wobei W = + V6/|a| ist, gilt r(w) = s(w). In den 


Intervallen — W < w < 0 und w > W ist demnach 
r(w) > s(w) , 
wihrend in den restlichen beiden Intervallen der w-Achse 


r(w) < s(w) 
gilt. Wegen 
e— Bar w/3 — |] (w= .0, G=.0) 


erhalten wir demnach 


q Ww feo) 
I(q) = I(W) + f s(w) w?™+! dw < ih e & w2m +1 dw + i oP © 402) diy (gi) 
W 0 0 

und weiterhin 


—4 =4 
I(—q) > J_s(w) wml dw aoe wert hes ag Ole 


also insgesamt 
a eres hal 
I(q) — I(—4) < | Be w?m+1 diy - ees —- | CP? 467 ei 
0 = 


Durch Auswertung der Integrale ergibt sich hieraus 


4 ep ~ (L(g) — 1(—4)) < € P(x) — 2 €P Iia(p) + 1; 
wobei ‘i 
a= Bq, p=pw=+iss (7.7) 
und 


hn(2) = Y (HI (7.8) 


ist. Da wir bei unserer Abschatzung das Integral I(q) durch das Integral lim I(q) ersetzt haben, 
Cha Mee} 

so gewinnen wir eine obere Schranke nicht nur fiir die Grenzamplitude Q*, sondern sogar fiir den 

Betrag einer jeden Amplitude Q,, n > 1. 


Das Ergebnis ist: Im Falle einer Riickstellkraft, die der Funktion 


of (qr early eae) 


entspricht, ist die Wurzel x =x, der Gleichung 
e* hy (x) = 2 eP hy(p) — 1 (7.9) 
eine obere Schranke fiir 6IQ,,| (n = 2, 3,...). Hierbei ist p durch (7.7) und h, durch 


(7.8) gegeben. 

Die Gleichung (7.9) kann z. B. mit Hilfe der Newtonschen Naherungsmethode behandelt werden, 
wobei man zweckmaBigerweise « = p als Ausgangswert wahlt. 

Der Fall der linearen Riickstellkraft entspricht dem Werte m = 0. Die Gleichung (7.9) hat 
dann die Form 


e* (x — 1) = 2eP (p—1) +1, (7.9') 


und dies laBt sich, da p stets positiv ist, durch die noch einfachere Gleichung 


e* (x — 1) = 2 peP (77) 
ersetzen. 

Das soeben behandelte Verfahren kann ohne weiteres auf den allgemeinen Fall ibertragen 
werden, wo f(q) ein ungerades Polynom mit positiven Koeffizienten ist. Man hat dann das In- 
tegral I(q) in cine Summe von Integralen zu zerlegen, deren jedes jeweils genau einen Summanden 
der Funktion f(q) im Integranden enthalt, und kann diese Integrale auf die soeben angegebene 
Weise abschatzen. Zu beachten bleibt lediglich, daB bei der Auswertung der Integrale verschiedene 
Potenzen von f auftreten und diese somit zum Teil in die zu (7.9) analoge Gleichung eingehen. 
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Wir betrachten hierzu den praktisch wichtigen Fall 


fM@=qtwv@. 


Da f(q) aus zwei Summanden besteht, ist I(q) in zwei Integrale aufzuspalten, und diese sind einzeln 
abzuschatzen. Als obere Schranke fiir 6 Q* erhalt man in diesem Falle die Wurzel x = x, der fol- 
genden Gleichung: 


e* h*(x) = 2 eP h*(p) + 62 + 6 pe; (7.10) 


h*(z) = B? (2 —1) +2 (2&—32+4 6z—6), 


und p ist durch (7.7) gegeben. Diese Gleichung lat sich ahnlich wie (7.9) lésen. 

Bisher haben wir ausschlieBlich echte Riickstellkrafte (vgl. (2.3)) betrachtet, da diese vor allem 
von praktischer Bedeutung sind. Wir zeigen nun, daB man auch im Falle eines beliebigen ungeraden 
Polynoms f(q) in ahnlicher Weise eine obere Schranke fiir Q* bestimmen kann. 

Man denke sich I(q) in eine Summe von Integralen zerlegt, deren jedes genau einem Glied des 
Polynoms f(q) entspricht. Es sei I+(q) bzw. I-(q) die Summe aller Glieder dieser Zerlegung, die 
(fir g > 0) einen positiven bzw. einen negativen Wert haben; offenbar sind das die Integrale, die 
einem Glied von f(q) mit positivem bzw. mit negativem Koeffizienten entsprechen. Man hat dann 


I(q) = '(q) + -(q)- 


Entsprechend zerlegt man I(— q). Es sei J+(—q) die Summe aller Integrale der Zerlegung von 
I(—gq) mit positivem Wert; offenbar sind dies nun diejenigen, die Potenzen von f(q) mit positiven 
Koeffizienten entsprechen. Es sei J~(—gq) die Summe der iibrigen Integrale der Zerlegung; diese 
besitzen einen negativen Wert. Es ist 


IN ict) ag) 
Man bestimme nun auf die vorstehend angegebene Weise Schranken A, B, C und D derart, daB 


Pq<A, J(-9g>B, —Fq@>e¢, —J(-g<D (7.11) 
gilt. So ergibt sich 


hierbei ist 


I(q) — I(—q) <x A— B— C+D. (7.12) 
Wiederum findet man die gesuchte obere Schranke fiir Q*, indem man die rechte Seite dieser Un- 
gleichung gleich null setzt. 

Der praktisch interessierende Fall 
aD Ae 
ist von dieser Art; f(q) besitzt auBer q = 0 noch die weiteren Nullstellen ¢g = — u-?. Dem Gliede 
p2 q® entsprechen die Bestandteile I+(q) und J*+(—gq) der obigen Zerlegung; dem Gliede —gq ent- 
sprechen die Bestandteile I~(q) und J~(—q). Man findet 
q as q 


— Fg) = | swede > [waw + | sw) wd > a= C 


und 
fan Bey —4q i 
— J~(—4q) = | s(e) w dw < | w dw + ee [ we au Gere © (3 ge c)=D, 
0 0 =¥ 
wobei 
V3 
Vane lees 
Tal 


ist. Indem man die dem Glied py? q® entsprechenden Integrale in der obenstehend angegebenen 
Weise abschatzt, erhalt man insgesamt eine obere Schranke fiir 8 Q* als Wurzel der folgenden Glei- 
chung in x: 


eo v) — ego (xt 24) = 260 lp) + 6 (7.18) 
mit ; 
ES pies = ans 22 = 
A= she g(z) = 2— 327 + 62—6 


und p gemaB (7.7). 
28 
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8. Verhalten des Schwingers fiir kleine Ausschlige. Bei kleinen Ausschlagen kann man Nahe- 
rungsaussagen iiber die Amplituden der Schwingungen dadurch erhalten, daB man eine endliche 
Anzahl von Gliedern der Taylorentwicklung der Funktion I(qg) betrachtet. Es ist 


sto) =14+ Bw + Fut +(F—e)ut + (Sp — c) tt fo (8.1) 


Durch Multiplikation dieser Reihe mit der Entwicklung der Funktion f(w) an der Stelle w =0 
und anschlieBende Integration ergibt sich die gewiinschte Potenzreihendarstellung fiir I (q). 
So findet man im Falle der Funktion 


fg=P"t* (m=0,1,...) 
[vgl. (6.2)] die Entwicklung 


acai Ba, (B2)a | BGP, 
ig =o"H (sara tim atime + Imt5 7 ). hae 


Wir bezeichnen die Summe der in (8.2) explizit angegebenen Glieder mit N(q). Wegen (3.4) ist die 
einem Anfangsausschlag q = Q, > 0 entsprechende nachste Amplitude, Q,, naherungsweise durch 


die Wurzel der Gleichung 
N(q) — N(Q,) = 9 (8.3) 


gegeben. Diese Gleichung lat sich am einfachsten mit einem der tblichen Naherungsverfahren, 
etwa dem Newtonschen, behandeln. Dabei wahlt man als erste Naherung fiir die gesuchte Wurzel 


zweckmaBigerweise den Wert q, = —Q,. Als zweite Naherung fiir Q, ergibt dann das Newtonsche 
Verfahren Ay Ron 
B/(2 m + 3) + (63/6 — c) Q3/(2 m + 5 

= f). =D GE 8.4 

=~ 20 TRA TP UR GIO” “7 


also speziell im Falle der linearen Riickstellkraft (m = 0) 


B/3 — (63/6 — c) Q2/5 
= — O20 + PER PH —)G * 


Die Naherung (8.4) gilt umso genauer, je kleiner Q, ist. Fiir nicht zu groBe Werte von |a| ist das 
Korrekturglied negativ. Dann wird |Q,| > Q,; man befindet sich in dem Gebiete der Energiezufuhr. 
Ist |x| hinreichend gro, so kann der Zahler des Korrekturgliedes selbst fiir kleine Werte von Q, 
einen negativen Wert annehmen. Dann ist |Q,| < Q,. In der Tat entspricht groBen Werten von |a| 
ein kleines Gebiet der Energiezufuhr, und (8.4) gestattet dann sinnvolle Aussagen tiber Q, auch bei 
Wahl von Werten Q,, die zwar relativ klein sind, aber trotzdem bereits im Gebiet des Energie- 
entzuges liegen. 


9. Die Grenzamplitude im Falle groBer Werte von x. Die soeben betrachtete Taylorentwicklung 
laBt sich auch zur Gewinnung von Naherungsaussagen iiber die Grenzamplitude Q* heranziehen. 
Damit man bei Beriicksichtigung einer geringen Anzahl von Gliedern dieser Entwicklung noch 
praktisch verwertbare Naherungen erhalt, mu8 Q* hinreichend klein sein; dies trifft offenbar fiir 
groBe Werte von |q| zu. 

Wir betrachten Schwingungen, die der Funktion 


fq@=¢rt* (m = 0,1,...) 
[vgl. (6.2)] entsprechen. Dann ist 


q 
Hg) 19) | sleds tdae a2 gi hte Ee 


= 
Die Grenzamplitude Q* ist naherungsweise gleich der Wurzel der Gleichung, die man erhalt, wenn 
man eine endliche Anzahl von Gliedern der Entwicklung (9.1) beriicksichtigt und die sich so er- 
gebende Funktion gleich null setzt. So erhalt man bei Verwendung der ersten beiden der genannten 
Glieder eine quadratische Gleichung, deren Lisung fiir gegeniiber « kleine Werte von f die Form 


se |) 2 ee) a 
Q Ly lL (2 m + 3) (m od 0, I, oe 5) (9.2) 
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hat. Wie man sieht, nimmt Q* mit wachsendem m monoton ab; dies ist vom Standpunkte der 
Mechanik aus einleuchtend. Ferner folgt aus (9.2) 


lim Qr = + 18 


(vgl. Abschnitt 4). 


10. Beispiele. Abb. 1 zeigt das »Amplitudendiagramm“ einer durch die Differentialgleichung (2.1) 
beschriebenen Schwingung mit linearer Riickstellkraft fiir den wesentlichen Parameter 
u=a/6 = 1. Dieses Diagramm wurde durch numerische Auswertung des Integrals I(q) gewonnen. 


Auf Grund der Amplitudenbeziehung 
Iq) 


(3.4) ist klar, wie man aus dieser gra- 
phischen Darstellung zu vorgegebener 
Anfangsamplitude Q, die Betrage der 
Amplituden Q,, Q3,...entnehmen kann. 

Ist zB. 8 =2 und wird Q,=1/2 7 
gewahlt, so findet man, wie in Abbil- 
dung 1 angedeutet ist, von BQ =1 
ausgehend nacheinander 


Q, & — 51,65 w — 0,83, 
ion anti 2 


usw. Wie man aus Abb. 1 ersieht, hat 
die Grenzamplitude Q* im vorliegen- 
den Falle etwa den Wert 2,2/8 = 1,1. 
Die Ungleichung (4.3) liefert die untere 
Schranke 0,86... Als obere Schranke 
erhalt man aus (7.9’’) den Wert 1,62 1 
und aus (7.9’) den (stets genaueren) 
Wert 1,45. Fir den Fall der linearen 
Rickstellkraft hat die Naherungs- 
formel (9.2) die Gestalt 

Q* & aes (9.2’) 


|2| 


und liefert in unserem Fall den Wert he : 
At, declan etwa 1%, eaten ge- Abb. 1. ee gee luk BOR gemaB (2.1) 
nauen Wert abweicht. 

Abb. 2 zeigt die Funktion 6 Q* in Abhangigkeit von dem wesentlichen Parameter u und zwar 
fiir den Fall f(q) = q (lineare Riickstellkraft; ausgezogene Kurve) und fiir den Fall f(q) = q? (ge- 
strichelte Kurve). Im erstgenannten Falle betragt die Abweichung der durch (9.2') gegebenen 
Naherungswerte von den exakten Werten bei dem in Abb. 2 wiedergegebenen Bereich etwa 1% 
bis 3%. Da die Naherungswerte auf einer Hyperbel liegen, so sind die exakten Werte, von diesem 
Fehler abgesehen, ebenfalls proportional zu 1/u. Auch fiir den Fall f(q) = q? ist die Genauigkeit der 
Naherungsformel (9.2) befriedigend. 


11. Behandlung der Differentialgleichung (1.3) mit Hilfe des Isoklinenverfahrens in der Phasen- 
ebene. Um die Differentialgleichung (1.3) mit Hilfe des Isoklinenverfahrens behandeln zu kénnen, 
ist es zweckmaBig, q = v, also q = (dv/dq) v zu setzen, also die Koordinaten q und v der Phasen- 
ebene zu verwenden. So erhalt man fiir negative Werte von v 


d Las 
fo +E (Lag) + Ff =0 (11.1) 
und fiir positive Werte von v 
d m, 2 = e : a 
Bo+ Fe + 5 (L@)—f@)=0 mit G=—a, (11.2) 


wobei f,, und f, wie in Abschnitt 2 definiert sind. Die Isoklinenbilder dieser beiden Differential- 
gleichungen lassen sich zu einem einheitlichen Isoklinenbild zusammensetzen, indem man die obere 


28" 
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(q, v)-Halbebene mit dem zu (11.2) gehérenden Isoklinenbild und die untere (q, v)-Halbebene mit 
dem zu (11.1) gehérenden Isoklinenbild zu einer Ebene zusammenfiigt. 


PQ" 


4 + ~ 


0 if 2 
Abb, 2. Dimensionslose Grenzamplituden £Q* in Abhangigkei ii 
é s ; gigkeit von «/f fiir 
Schwinger gemaB (2.1) mit f(q) = q (ausgezogen) und f(q) = q° (gestrichelt), 


Abb. 3 zeigt ein solches kombinier- 
tes Isoklinenbild fiir den Fall der line- 
aren Riickstellkraft und die Parameter- 
werte « =p = 1. Die Isoklinen 
haben in diesem Fall fiir (11.1) die 
Darstellung 
v=(P—1)7 (b+ P+ 4q¢—)), 
wobei dv/dq = k = konst. gesetzt ist; 
hierbei sind jeweils nur diejenigen 
Kurvenzweige zu beriicksichtigen, fiir 
die v < 0 ausfallt. Eine ahnliche Dar- 
stellung ergibt sich fiir die Isoklinen 
der Differentialgleichung (11.2). 

In Abb. 3 ist die Grenzschwingung 
(gestrichelt) eingezeichnet und weiter- 
hin eine bei g = Q, = 0,25 beginnende 
Schwingung, deren Amplituden dem 
Betrage nach monoton der Grenzampli- 
tude zustreben. 

Wegen v = dq/dt ist 


(11.3) 


und man kann, indem man die Werte 
des Integranden dieses Integrals aus 
dem Isoklinenbilde entnimmt, durch 


tL tL fecprtrc oder numerische Inte- 
4 f gration die Funktion q = q(t) fiir die 


Schwingung 


jeweils interessierende 
naherungsweise ermitteln. 


v 
a 


Bs 


H l A 


—a 


=! 


-3 -2 =7, 0 


a 3 
i 


Abb. 3, Isoklinenbild, Integralkurve (fiir Q, = 0,25) und Crenzzykel (gestrichelt) fiir Schwinger gemaB (1.3) mit f@= undo =pP=y=1 
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12. Eine Naherungslésung der Differentialgleichung (1.2) fiir kleine Werte von fb und lineare 
Rickstellkraft. Wir betrachten die Differentialgleichung (1.2) wiederum fiir eine aus einer Ruhe- 
lage q = Q, > 0 beginnende Schwingung. Dann ist bis zum nichsten Umkehrpunkt q = Q, < 0 
z lie ' 
stets q<. 0. Mitd = oi: 6 nimmt (1.3) fiir diesen Abschnitt der Schwingung im Falle linearer Riick- 
stellkraft die Form 
g+6(l—o8@) @+22q=—0 (12.1) 
an. Wir untersuchen die Lésung dieser Differentialgleichung fiir kleine Werte von 6. 
Fir 6 = 0 lautet die allgemeine Lésung 


q(t) = a sin (xt + b) (12.2) 
und die der obigen Anfangsbedingung entsprechende Sonderliésung 
q(t) =Q, cosxt. (12.2") 


Fir kleine Werte von |6| wird sich die Lésung von (12.1) nur ,,wenig von (12.2) unterscheiden. 
Demnach ist zu erwarten, daB man fiir kleine Werte von |d| mit Hilfe der ,,Methode der langsam 
variablen Amplitude und Phase“ (d. h. unter der Annahme, daB die in (12.2) auftretenden GréBen a 
und 6 zwar zeitabhangig sind, sich aber nur ,,langsam“ andern) eine brauchbare Naherungslésung 
erhalt. 

Wir bemerken jedoch sogleich, daB die itibliche Form dieser Methode?, wie man sie z. B. bei 
der van der Polschen Differentialgleichung verwenden kann, hier nicht zum Ziele fiihrt. Man kommt 
namlich bei diesem Grade der Approximation wieder zu der Ausgangsannahme, da a und 6 kon- 
stant sind, also zu (12.2) zuriick. Dies ist dadurch bedingt, da die Funktionen q und ¢ im Dampfungs- 
glied von (12.1) quadratisch vorkommen. Wir werden die bei Anwendung der tiblichen Methode 
durchzufiithrende Mittelwertbildung durch eine andere geeignete Integration ersetzen; diese schlieBt 
eine geringere Vernachlassigung in sich und ergibt eine Lésung, die (12.1) bis auf Glieder in 6? und 
héherer Ordnung befriedigt. 

Differentiation von (12.2) nach t liefert 


P—axcosp mit pox b. (12.3) 
Wir haben eine Lésung von der Form 


q = a(t) sin p (12.4 


anzugeben und die nun zeitabhangig zu denkenden GréBen a und 6 zu bestimmen. Durch Differen- 
tiation folgt aus (12.4) 


q = asin p +a(x +6) cosp. (12.5) 
Indem wir den durch (12.3) gegebenen Naherungsausdruck in (12.5) einfiihren, ergibt sich 
asin p +abcosp=0. (12.6) 
Fassen wir in (12.3) nun die GréBen a und 6 als zeitabhangig auf, so finden wir durch Differentiation 
q = ax cos p—ax (x + b)sinp. (12.7) 
Setzen wir (12.2), (12.3) und (12.7) in (12.1) ein, so erhalten wir 
a cos p —ab sin p = — 6 a? x (1 — «7 a’ sin? p) cos® p. 
Unter Beriicksichtigung von (12.6) laBt sich diese Gleichung in das folgende System von Differen- 
tialgleichungen zerlegen: 


a = —6 a?x (1 —o? a? sin? p) cos* p , (12.8) 


b = 6 ax (1 — «x? a? sin? p) cos* p sin p 
Nach der iiblichen Methode der langsam variablen Amplitude und Phase ware nun der Mittelwert 


der beiden linksstehenden Groen iiber eine Periode zu bilden. Dies ergibe a fs 0 und b=0. Wir 
integrieren stattdessen die Differentialgleichungen (12.8); hierbei sehen wir die in den rechten 


1 N. W. McLachlan, Ordinary non-linear differential equations. Oxford 1950. 
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Seiten enthaltenen GréGen a und b als konstant an. So ergibt sich 


. 3 2 5 ; 
a(t) = a — 6a? (sin p — (1 + a? a2) = P L «2 a? ~ | (12.9) 


und 


5 
b(t) =b—da((l a? a?) FP 4. 92 4? sek (12.10) 


wobei die (rechts an erster Stelle stehenden) Integrationskonstanten wegen (12.2) mit a baw. b 
bezeichnet worden sind. Hieraus lassen sich nun leicht die Fourierentwicklungen von a(t) und b(t) 


il : i 1 Swit 
gewinnen; indem wir 6 wieder durch 3 P ersetzen und die Abkiirzung A = zo a* einfiihren, er- 
halten wir 


a(t) = a —(f a/8) (@ m2 A) sin p(t ed Sa) . (12.11) 


b(t) = b — (8 a/8) (a 2A) cod pi (Ut) ’) ; (12.12) 


Bemerkenswerterweise bestehen diese Entwicklungen aus nur endlich vielen Gliedern. 

Betrachten wir in (12.11) und (12.12) nur die langsam veranderlichen Bestandteile, also nur 
die sin p baw. cos p enthaltenden Glieder, so gewinnen wir die Naherung 

q(t) = (a — r(a) sin p) sin (% t + b — z(a) cos p) (12.33) 
mit ; 
It 
r(a) == fa (3—2 4), (a) = 7h a(1—2 4). 
Ly 2 

Setzen wir die der Lésung (12.2') entsprechenden konstanten Werte a = Q, und b = 77 in (12.13) 
ein, so ergibt sich die Sonderlésung 


q(t) ~ (Qi — r(Q,) cos x t) cos (x t + 2(Q,) sin x t) , (12.14) 


q(9) = Q —7(Q) » q0) = 0" 


Dies besagt, da den Anfangsbedingungen fiir (12.2) nur geringfiigig geanderte Anfangsbedingungen 
fiir (12.13) entsprechen, indem sich die Ausgangslage von Q, nach Q, — r(Q,) verschiebt. Weiterhin 
sieht man schlieBlich, daB r(a) fiir 


also fiir t = 0 speziell 


petal. (12.15) 


verschwindet. Dieser Wert stellt demnach eine Naherung fiir die Grenzamplitude dar [vgl. (12.14)]; 
er weicht allerdings von der durch (9.2') gegebenen Naherung (s. Abschn. 10) um rund 10% ab. 


13. Zusammenfassung. Die betrachtete ,,modifizierte van der Polsche Differentialgleichung“ 
@ — (sgn) & 1 —o8 @) P+ 2 f(@) = 0 (1.3) 


beschreibt das Verhalten einer besonderen Klasse von selbsterregten Schwingern mit einem Frei- 
heitsgrad. Die Gleichung (1.3) 1aBt sich in eine lineare Differentialgleichung erster Ordnung iiber- 
fihren; mit V = q? und y? = 2+? erhalt man im Falle negativer Geschwindigkeit 


Vi + B(L—o @)V + fq) =0; ere 


im Falle positiver Geschwindigkeit ergibt sich eine ahnliche Differentialgleichung [namlich (2.2’)]. 
Die Riickstellkraft des Systems ist zunachst beliebig. Vom physikalischen Standpunkt aus ist es 
sinnvoll, sich auf ,,echte“* Riickstellkrafte zu beschranken, d.h. auf Rickstellkrafte, die stets im 
Sinne einer Verringerung des jeweiligen Ausschlages wirken. Die Funktion f (q) wird deshalb den 
Bedingungen 


S@ 20 fir q>0, F@)S0 firg=0., (q) = 0 (2.3) 

unterworfen. . 
Die Lésung der Differentialgleichung (2.2) 1aBt sich, von trivialen Sonderfallen abgesehen, 
nicht durch eine endliche Anzahl bekannter Funktionen darstellen. Trotzdem kann man weit- 
gehenden Aufschlu8 iiber das Verhalten des Schwingers erhalten, insbesondere auf Grund der 
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»,Amplitudenbeziehung“ (3.4), mit deren Hilfe sich die Folge der Umkehrpunkte der Schwingung 
bestimmen laBt. 

Die in der vorliegenden Untersuchung hergeleiteten Ergebnisse beziehen sich auf eine Auswahl 
praktisch wichtiger Formen der Riickstellkraft. Wie man aus den Betrachtungen ersieht, gestatten 
die verwendeten Methoden die Herleitung ahnlicher Resultate auch fiir andere Riickstellkrafte. 
Wir nehmen stets « ~ 0 an; der Falla = 0, der von dem Fall x ~ 0 vom mechanischen Standpunkt 
aus grundsatzlich verschieden ist, wurde frither! schon ausfiihrlich untersucht. 

Fiir kleine Ausschlage q wird dem Schwinger Energie zugefiihrt, wihrend fiir hinreichend groBe 
Werte q Energie entzogen wird. In Abschnitt 5 wird bewiesen, da®, wenn f(q) ungerade ist, die 
Folge der Betrage der Maximalausschlage Q,, Q.,... (im Text kurz als ,,Amplituden“ bezeichnet) 
gegen einen endlichen Grenzwert Q* (,,Grenzamplitude“) strebt. Die genannte Folge wachst 
bzw. fallt monoton je nachdem, ob die Schwingung von einer Ruhelage q = Q, < Q* oder 
q =Q, > Q* aus beginnt. Ist f(q) nicht ungerade, so erhalt man statt des einen Grenzwertes Q* 
zwei Grenzwerte. 

Die GréBe der in der Differentialgleichung (2.2) auftretenden Parameter x, f und y ist von 
unterschiedlichem Einflu8 auf das Verhalten des Schwingers. Ist beispielsweise f(q) = q?"+}, 
m = 0, 1,..., so hangt das Verhaltnis der Betrage aufeinanderfolgender Amplituden nur von dem 
»,wesentlichen Parameter“ u = «/8 ab und wachst monoton mit wachsendem u (vgl. Abschnitt 6). 

In den Abschnitten 4 und 7 werden praktisch brauchbare untere und obere Schranken fiir die 
Grenzamplitude Q* angegeben. Zum Beispiel gilt die Abschatzung Q* > 7/3/|a|. Obere Schran- 
ken ergeben sich als Nullstellen gewisser Gleichungen [vgl. (7.9) und (7.10)]. Derartige Schranken 
lassen sich auch im Falle gewisser Riickstellkrafte gewinnen, fiir die (2.3) nicht erfillt ist (vgl. im 
zweiten Teil des Abschnittes 7). 

Fir kleine Ausschlige kann man Naherungsaussagen iiber die Amplituden erhalten, indem 
man eine endliche Anzahl von Gliedern der Taylorreihe der. Lésung von (2.2) beriicksichtigt (vgl. 
Abschnitt 8). Dieser Weg fiihrt auch zu (bei groBen Werten von |x| praktisch brauchbaren) Nahe- 
rungen fiir die Grenzamplitude Q* (vgl. Abschnitt 9). 

In Abschnitt 11 wird die betrachtete Differentialgleichung mit Hilfe des Isoklinenverfahrens 
in der Phasenebene behandelt. SchlieBlich wird in Abschnitt 12 eine Naherungslésung der Dif- 
ferentialgleichung (1.3) fiir den Fall kleiner Werte von # und lineare Riickstellkraft auf eine Weise 
hergeleitet, die eine Verfeinerung der iiblichen ,,Methode der langsam variablen Amplitude und 
Phase“ darstellt. 

(Eingegangen am 17. Januar 1957) 


Anschriften der Verfasser: 
Professor Dr.-Ing. K. Klotter, Stanford, Calif., Division of Engineering Mechanics, Stanford University 
Professor Dr. E. Kreyszig, Columbus, Ohio, Department of Mathematics, Ohio State University 


1 Siehe FuBnote 5 von Seite 389. 
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Piezoelektrische Gerate zur Bestimmung des Pitotdrucks, 
des statischen Drucks, des dynamischen Drucks und des Impulses 
einer pulsierenden Stromung 


Von F. Staab 


1. Einleitung. Man kann ohne Schwierigkeit, z. B. in Windkandlen, bei stationaren Str6mungen 
mit Hilfe eines Réhrchens (Pitotrohr), das der Strémung entgegengestellt wird und durch welches 
man den Druck der zur Ruhe gekommenen Gase erhalt, bei Unter- und Uberschallgeschwindig- 
keiten auf die Gasgeschwindigkeiten schlieSen. Diese Methode liefert bei pulsierenden Str6mungen 
wohl Mittelwerte, doch kann daraus kaum der wirkliche veranderliche Strémungsvorgang erhalten 
werden. 

Es wurde jedoch viel Miihe aufgewendet um bei pulsierenden Strémungen die durchflieBenden 
Gasmengen, etwa durch ein Rohr, mit Hilfe von Diisen und Blenden zu erhalten } 2,3,4 und da- 
durch auch die mittleren Geschwindigkeiten. Ebenso gibt es eine gréBere Zahl von MeBgeradten 
um den statischen Druck in einer stationdaren Strémung durch ein Hakenrohr, das seitliche Offnun- 
gen hat, oder bei pulsierenden Strémungen am Rohrrand durch piezoelektrische MeSelemente 
zu bestimmen ® ® 7:8, Durch diese MeBanordnungen erhalt man am Rande der Strémung mit 
ziemlicher Genauigkeit bis zu héheren Frequenzen den verdnderlichen statischen Druck. 

Man hat auch Prallplatten benutzt, um mit ihnen die RiickstoBkrafte z. B. eines Motors, wenn 
dessen Strémung gleichférmig verlauft zu bestimmen und diese Methode auch fiir Pulsationen 
angewandt, wobei jedoch nur ein Mittelwert des Impulses und des Schubes bestimmt wurde’. 

Bei den Arbeiten an pulsierenden Strahltriebwerken in Deutschland wahrend des Krieges ergab 
sich die dringende Notwendigkeit, deren rasch veranderliche Gasgeschwindigkeiten zu bestimmen, 
ebenso den statischen Druck im Strahl, da Messungen an den Triebwerken selbst oft von heftigen 
Erschiitterungen gestért waren. Ebenso war es notwendig, den veranderlichen Impuls im Strahl 
zu erhalten, da Schubmessungen bei rascheren Verbrennungen etwa kaum médglich waren. 
Aus diesen Griinden wurde von 1941 an eine Zahl von MeBgeraten entwickelt, wobei der piezo- 
elektrische Effekt beniitzt wurde, um die genannten GréBen zu erhalten. Diese MeBgerate haben 
sich tiber Jahre bewahrt. 

Schadliche instationare Strémungszustande treten auch bei den modernen Strahltriebwerken auf. 
Bei der Nachverbrennung in Strahlturbinen z. B., ebenso bei Raketen und Staustrahltriebwerken. 
Die Erfassung der instationéren turbulenten Ausbreitungsvorgange im Strahl ist dazu wegen der 
damit vorhandenen Larmaussendung, die verringert werden muB, sehr wichtig. Ebenso sind MeB- 
gerate dieser Art fiir gewisse strémungstechnische Untersuchungen von gréBter Bedeutung, und 
es ist von groBem Wert, bequeme und sicher arbeitende Untersuchungsanordnungen zu besitzen. 

Es ist auch nur schwer méglich, mit Hilfe der Schwebeteilchenkinomatographie leuchtende 
Teilchen festzuhalten, Ionenlaufzeitmessungen durchzufiihren, ebenso die Schlierenkinomato- 
graphie anzuwenden und StoBwellen, weiter Ultraschallwellen aufzunehmen bei den rasch ver- 


anderlichen Vorgangen, die sich auch iiber den Querschnitt des Strahls, z. B. von Strahlantrieben, 
rasch verdndern 20, 11, 12, 13, 


1 Q. Lutz, Ing.-Arch. 3 (1932) S. 138 und S. 432. 

> F. Schultz-Grunow, Forschung Ing.-Wesen 12 (1941) S. 117. 

° A. K. Oppenheim und E. G. Chilton, Transactions ASME (1955) S. 231. 

4 W. A. Wildhack, Science 120 (1954) S. 191. 

° S. Meurer, Forschung Ing.-Wesen 11 (1940) S. 237. 

® K. Staiger, Motortechn. Z. 5 (1944) S. 87. 

* W. Gohlke, Einfiihrung in die piezoelektrische MeBtechnik, Leipzig 1954. Mit einer umfassenden Schrift- 
tumsangabe, welche besonders auch die wichtigen Arbeiten des Verfassers enthalt, die ehemals im Institut fiir 


Pg EEE der Luftfahrtforschungsanstalt Braunschweig (Leiter Prof. Dr. E. Schmidt) durchgefiihrt 
wurden. 


8 T. Wrathall, Instruments 26 (1953) S. 50. 

° S. Hasinger, Jb. 1940 d. Dtsch. Luftfahrtforschung, S. II 304. 

‘© H. Schardin, Ergeb. exakt. Naturw. 20 (1942) S. 303. 

me High Speed Aerodynamics and Jet Propulsion IX (1954) Princeton University Press: Physical measure- 
ments in Gasdynamics and Combustion. R. W. Ladenburg, B. Lewis, R. N. Pease, H. S. Taylor. 

12 G. Dubois und R. Kling, C. R. Acad. Sci. Paris 228 (1949) S. 363. 


18 D.G. Marlow, C. R. Nisewanger und W. M. Cady, J. appl. Phys. 20 (1949) S. 771. 
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Da man dazu Ionen, deren Menge sich bei pulsierenden Vorgangen rasch verandert, in den 
Strahlen von Strahlantrieben antrifft, ist es dazu kaum méglich, mit Hilfe der Koronaentladung 
Gasgeschwindigkeitsbestimmungen durchzufiihren!. Schwierig ist weiter die Anwendung von 
Hitzdrahtanemometern?. Wenn man auch bei Luftstrahlmodellen geringen Querschnitts brauch- 
bare Ergebnisse mit Hilfe von Photo-Multipliern gewonnen hat, so ist diese Methode doch be- 
schrankt bei den gedachten Anwendungen. 

Es ist dazu bequem, Verstarker und Oszillographen zu beniitzen, die an und fiir sich schon zur 
Bestimmung des statischen Druckes am Rohrrand vorhanden sind. 


2. Prinzip der Mefanordnungen. Man kann den Pitotdruck in einer pulsierenden Strémung 
bestimmen, bei der sich die Gasgeschwindigkeiten und der statische Druck rasch andern, indem 
man die Drucke iiber ein kleines kolbenartiges Gebilde auf Quarze wirksam sein laBt. Auf deren 
Oberflachen entstehen sich andernde Elektrizitétsmengen. AuBerlich sieht die MeSeinrichtung 
ahnlich aus wie MeBgerate dieser Art fiir stationare Strémungen. Ein Verstarker erhéht die Span- 
nungsainderungen, auch eines zumeist parallel geschaltenen Kondensators, die mit Hilfe von 
Schleifen- oder Kathodenstrahloszillographen angezeigt werden. Die hohen Eigenfrequenzen der 
Anordnungen erlauben noch betrachtliche Frequenzen aufzunehmen. , 
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Abb. 1. Gerat zur Bestimmung des Pitotdruckverlaufes in einer pulsierenden Strémung. Abb. 2. Hohlzylinderquarz, 
1 Hohlzylinderdoppelquarz, 2 Réhrenfeder, 3 Endteller, 5 Vorspannschraubenmutter, der aus zwei Halften besteht. 


6 und 7 Kihlwasserleitungen, 8 Kupfermembran, 9 Gehause, 10 Vorderteil, 11 Metallfeder, 
12 Leitung zum Verstirker. 


Zur Bestimmung des veranderlichen statischen Druckes im Strahl beniitzt man ein Gerat, das 
ahnlich wie ein Hakenrohr aussieht, mit scharfer oder abgerundeter Spitze, hinter der in einem 
gewissen Abstand, so daf die Stérungen, welche von der Spitze ausgehen den veranderlichen stati- 
schen Druck nicht mehr merklich beeinflussen, Schlitze angebracht sind. Diese miissen so breit 
sein, daB der statische Druck zu dem kolbenartigen Gebilde, welches die Krafte aufnimmt, dringen 
kann. Auf ahnliche Weise erhalt man weiter den veranderlichen Impuls einer Strémung. 


3. Geriite zur Bestimmung des Pitotdrucks. Der Aufbau eines Druckelementes der genannten 
Art kann folgendermaBen geschehen*. In Abb. 1 stellt 1 einen Hohlzylinderquarz dar. Dieser ist 
aus zwei gleich groBen Halften gebildet, die so zusammengelegt sind, daB die elektrischen und die 
neutralen Achsen in entgegengesetzte, die optischen Achsen in die gleichen Richtungen zeigen, 
wie dies in Abb. 2 angedeutet ist. Bei Druck in Richtung der Hohlzylinderachse entstehen auf 
den duSeren Flachen negative Ladungen, wahrend sich innen positive bilden. Die auSeren und in- 
neren Flachen der Quarze sind mit einer Metallschicht versehen und die inneren geerdet. An den 
Enden der du eren Zylinderflachen der Quarze sind schmale Bereiche, die als Isolationsflachen 
dienen >. 


1 W. Fucks, Z. Physik 137 (1954) S. 49. 
2 U. Jacob, Arch. techn. Messen 5 (1954) S. lu. 25. Vgl. auch S. 219—276 des in FuBnote Il von Seite 404 


genannten Buches. 


3 A, Powell, Proc. Phys. Soc. (B) 66 (1953) S. 1039. 
4 F. Staab, Gerate zur Bestimmung des Pitotdrucks einer pulsierenden Stromung, FKFS-Bericht 1942. 


5 Siehe FuBnote 5 von Seite 404. 
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Durch den Doppelquarz geht eine Réhrenfeder aus méglichst leichtem Werkstoff 2 hindurch 
(Abb. 1), welche mit ihrem Endteller 3 die Quarze gegen ein glattes Widerlager 4 mit Hilfe einer 
Schraubenmutter 5 driickt. Auf diese Weise wird eine gréBere Vorspannung erzeugt. In der 
Réhrenfeder wird durch das Réhrchen 6 Kihlwasser nach 3 gefiihrt. Dieses wird zwischen diesem 
Réhrchen und der Innenwand der Réhrenfeder bei 7 abgefiihrt. 

Zwischen dem Endteller 3 und den Doppelquarzen 1 ist eine ausgegliihte Kupfermembran 8 
von 0,1 mm oder einer geringeren Dicke eingepreBt, die von dem vorderen Teil 10 des Gehauses 
auf den Gehauseeinsprung 9 angepreBt wird. Auf sie wirkt eine verhaltnismaBig kleine Kraft, 
so daB rasche Temperaturschwankungen von ihr kaum die Anzeigen des Druckelementes beein- 
flussen. Der Vorderteil 10 kann halbkugelférmig oder kegelig ausgebildet sein und besitzt eine 
Bohrung, z. B. entsprechend den Massen, die beim Prandtlschen Staurohr iiblich sind. 

Auf den auBeren Flachen des Quarzes liegt noch eine aufgeschlitzte Metallfeder 11. Eine auf 
der Feder elastisch sitzende oder mit ihr verbundene Leitung 12 fiihrt zum Verstarker. Bei einer 
Anderung des Druckes und der mechanischen Belastung des Kristalls andert sich der Anoden- 
strom, wobei dessen Anderung ein MaB fiir die entstehende Ladungsanderung ist. Der riickwartige 
Gehauseteil wird durch eine entsprechend geformte Schraube 13 abgeschlossen. In den Gehause- 
raum mit den Quarzen kénnen keine Verbrennungsgase eindringen. 

Die bei einer statischen Eichung, etwa mit Druckluft, sich ergebende Eichkurve, also die Ab- 
hangigkeit des Oszillographenausschlages von der Gréfe des Druckes im Innenraum bei 10, ist 
linear zu erhalten. Die Empfindlichkeit war selbst bei langerem Kabel zwischen Quarzelement 
und Verstarker so groB, daB etwa fiir eine halbe Atmosphare Uberdruck im Innenraum bei 10 
einige cm Oszillographenausschlag erhalten wurden. Im iibrigen wurden iibliche Netzanschluf- 
verstarker verwendet, welche fiir die Bestimmung des Verlaufs und der Gréfe des statischen Drucks 
am Rohrrand mit Hilfe iiblicher Druckgeber beniitzt wurden}. 
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Abb. 3, Gerat zur Bestimmung des Pitotdruckverlaufs. Links: mit ebener Vorderflache. Abb. 4. Verbindung von zwei Hohlzylinder- 
Rechts: vorderer Teil mit Durchbohrung und scharfer Vorderkante. doppelquarzen, 


Es wird hier der transversale piezoelektrische Effekt verwendet. Wahrend beim longitudinalen 
piezoelektrischen Effekt die frei werdenden Elektrizitatsmengen von den Abmessungen des Kristall- 
stiickes unabhangig sind, ist dies beim transversalem Effekt nicht der Fall. Lat man den Pitot- 
druck G, in Richtung der y-Achse, der neutralen Achse, wirken, und sind J und d die Abmessungen 
eines Quarzstabes in Richtung der y- und x-Achse, also der elektrischen oder polaren Achse?, dann 
ist die frei werdende Elektrizitatsmenge 


Q, =—dy aa G, (1) 


wobei d,, der entsprechende piezoelektrische Modul ist. Man kann also durch eine Verkleinerung 
der Quarzdicke und eine Verlangerung des Quarzes die Empfindlichkeit erhéhen. 
Da die Elektrizitatsmenge, welche beim longitudinalen Effekt frei wird — G, ist wieder der 
wirksame Pitotdruck — 
Q, = dy 6, (2) 


ist, kann man dies hierbei nicht erreichen, und es lassen sich beim transversalen Effekt gréBere 
Empfindlichkeiten gewinnen. 

Druckelemente nach Abb. 1 wurden mit 10 mm Durchmesser verwendet. Doch kann dieser 
noch etwas verkleinert werden. Dieser Durchmesser ist verhaltnismaBig gering gegeniiber den 
gréBeren des Strahls der modernen Strahlantriebe. 

In Abb. 3 sind noch zwei andere Ausfiihrungen des vorderen Gerateteils zur Bestimmung des 
Pitotdruckes angedeutet, unten eine auch sonst iibliche Anordnung mit scharfer Vorderkante, 
die ebenfalls in das Gehaduse eingeschraubt werden kann. Es ist wertvoll, sich auch hier den tiblichen 


1 Siehe FuBnote 6 von Seite 404. 


; L. Bergmann, Der Ultraschall und seine Anwendung in Wissenschaft und Technik, 6. Aufl. I. Teil 
Kapitel 2. Stuttgart 1954. i 
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Formen anzupassen, um die Gréfe der Fehler abschatzen zu kinnen! 2 3, Danebenist eine fiir manche 
Zwecke interessierende Anordnung, bei welcher der kleinere Endteller 3 in Abb. 1 vergroBert ist. 
So wird sein duferer Durchmesser gleich dem AuSendurchmesser des Gerategehauses. D6ch kénnen 
hier auch Abrundungen usw. angebracht sein. Der Spalt zwischen diesem vorderen Teil und dem 
Gehause ist so breit, daB dazwischen der verdanderliche statische Druck der Strémung geniigend 
rasch eindringen kann. Dadurch erhalt man eine gréBere Empfindlichkeit der MeBSeinrichtung 
fiir den dynamischen Druck, doch auch unter Umstianden Nachteile. 

Um die Empfindlichkeit zu erhéhen, wurden zwei oder auch mehr Hohlzylinderquarze, die aus 
je zwei Halften bestehen von etwa 6 mm auSerem, 3,5 mm innerem Durchmesser und etwa je3cm 
Lange iiber einer Réhrenfeder zusammengeschoben und eine Vorspannung ausgeiibt. Dabei liegt 
zwischen ihnen ein Metallring, um die Kraft des Endtellers auch auf den zweiten Quarz sicher 
und gleichférmig zu iibertragen. Elektrisch werden die Quarze so gekoppelt, daB die sehr diinnen 
Metallagen der auBeren Quarzflichen durch einen elastischen und aufgeschnittenen Metallring 
verbunden werden. Die inneren leitenden Flachen sind wie beschrieben durch die Réhrenfeder, 


auf der diese anliegen, geerdet. Abb. 4 bringt eine Skizze. Die iibrige Anordnung ist ganz ahnlich 
wie bei den schon beschriebenen Geriten. 


4, Auswertung der Versuchsergebnisse. Man kann die Ergebnisse, die bei stationdren Stro- 
mungen gewonnen wurden, auch fiir diese Gerate zur Bestimmung von Pitotdruck verwenden, 
also die Hinfliisse der verschiedenen Formen fiir die unterschiedlichen Gasgeschwindigkeiten 
beniitzen, ebenso etwa die Abhangigkeit der Anzeige vom Winkel der Strémungsrichtung zur 
Gerateachse. Weiter die bekannten Auswertemethoden um die Gasgeschwindigkeiten zu erhalten 
fiir die einzelnen Zeitelemente anwenden * °: 6: 7, 

Bei Unterschallgeschwindigkeiten, bei der die Strémungsgeschwindigkeiten kleiner als die 
Laval- oder kritischen Geschwindigkeiten sind, mift man beim Aufstau mit einem Gerat dieser Art 
den der kinetischen Energie entsprechenden Druck (Kesseldruck). Bei Uberschallgeschwindigkeiten 
tritt vor dem Pitotrohr ein VerdichtungsstoB auf, und man erhdlt einen niedereren Druck. 

Um aus den Pitotdrucken etwa die Gasgeschwindigkeiten erhalten zu kénnen, sind noch die 
Molekulargewichte der Gase, die Temperaturen und die statischen Drucke in der ungestérten 
Strémung nétig. Diese kénnen mit den in folgenden beschriebenen Geraten bestimmt werden. 
Im allgemeinen kann man bei den verschiedenen Strahlantrieben das Molekulargewicht ziemlich 
genau abschatzen, auch Rauchgasanalysen durchfiihren. 

Bei pulsierenden Strahltriebwerken wuBte man z. B. aus Flammenaufnahmen, da beim Aus- 
stoB der Gasmassen in jeder Periode zuerst Luft ziemlich gewéhnlicher Temperatur ausgeschoben 
wird, dann erst leuchtende Gase. Tragheitslose Bestimmungen der veranderlichen Temperaturen 
sind hier méglich mit Hilfe von Photozellen, indem man etwa wenig Kalium dem Kraftstoff bei- 
mischt. Dabei geht man davon aus, da thermodynamisches Gleichgewicht in den leuchtenden 
Gasen vorhanden ist’. Doch kann dies angenommen werden, da selbst bei raschen Ausstr6mungs- 
vorgangen, z. B. aus Raketenbrennkammern, jeweils der Gieichgewichtszustand praktisch er- 
reicht wird ® 1°, 

An allen Strahlantrieben treten instationaire Zustande auf, bei Strahlturbinen besonders in der 
Nachyerbrennung, bei Raketen und Staustrahltriebwerken. Doch 1aBt sich eine Bestimmung der 
Pitot- und statischen Drucke auch bei weitgehend stationarem Verlauf der Gasgeschwindigkeiten 
iiber die Querschnitte so durchfiihren, da man vom Rand des Strahls aus zur Achse die MeBgerate 
hereinfiihrt. Eine nicht zu hohe Einfiihrungsgeschwindigkeit stért wenig, da die Strahlgeschwindig- 
keiten im Vergleich dazu wesentlich héher sind. Dazu kénnen diese Stéreinwirkungen auch rech- 
nerisch erfaBt werden. 

Die Eigenschwingungszahlen von Hohlzylinderanordnungen der hier verwendeten Art sind 
verhaltnismaBig hoch. Versuchsergebnisse dariiber, die mit Hilfe von Ausschwingkurven erhalten 


W. Wien und F. Harms, Handbuch der Experimentalphysik Bd. IV, 1. Teil S. 489. Leipzig 1931. 
0. Walchner, Jb. 1938 d. Dtsch. Luftfahrtforschung, 5S. I 578. 
D. J. Kettle, J. Roy. Aeron. Soc. 58 (1954) S. 835. os 
W. Wien und F. Harms, Handbuch der Experimentalphysik Bd. IV, i. Teil S. 376. Leipzig 1931. 
K. Oswatitsch, Gasdynamik, Wien 1952. } 
Hiitte Bd. I, 28. Aufl.: A. Betz, Mechanik elastischer Fliissigkeiten S. 825. Berlin 1955. Net 
7? E. Schmidt, Einfihrung in die Technische Thermodynamik. 6. Aufl. Kap. XIII. Berlin/Gottingen/ 
Heidelberg 1956. 
8 H. Graff, Luftfahrtforschung 18 (1941) S. 8. 
9 QO. Lutz, Z. Flugwiss. 3 (1955) S. 151. ; ; 
10 High Speed Aerodynamics and Jet Propulsion IT (1955) Princeton University Press: Combustion Pro- 
cesses, B. Lewis, R. N. Pease, H. S. Taylor, S. 23. 
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wurden, liegen nicht mehr vor. Doch kann man nach den vorliegenden Erfahrungen schlieBen?, 
daB Vorginge von vielen Hertz noch gut erfabt werden kénnen. ’ i 

Es ist ein geringer Zeitunterschied vorhanden zwischen dem Eintreffen einer Druckanderung 
an der Spitze des MeBgerates zur Bestimmung des Pitotdruckes gegeniiber dem Ankommen dieser 
Druckanderung beim Endteller der Rohrenfeder. Doch ist dieser Zeitunterschied, entsprechend 
der Raschheit einer Druckausbreitung, sehr gering, wenn der Weg von der Spitze zum Endteller 
nicht zu lang ist. Man wird diesen auch deshalb nicht zu groB wahlen, um Reibungs- und Tur- 
bulenzeinfliisse, ebenso die Einwirkungen von Schwingungen, médglichst klein zu halten. 

Eine statische Eichung mit Hilfe eines Gases, das mit verschiedenem Druck auf das Quarzgerat 
gegeben wird, ist zumeist méglich, da dessen Eigenschwingungszahlen verhaltnismaBig hoch liegen 
gegeniiber dem aufzunehmenden Vorgang. AuBerdem sind Ladungsverluste durch einen hohen 
Isolationswiderstand des Eingangskreises im Verstarker klein, so daB geniigend groBe Zeitkonstan- 
ten erhalten werden. Bei den Mefgeraten mit ebener Vorderflache ist auch naherungsweise eine 
Eichung mit Gewichten moglich, Doch kann eine solche im Windkanal wertvoll sein. Ks ist dazu 
méglich, eine Kichung mit Hilfe von Flammenaufnahmen durchzufihren auch bei veranderlichen 
Zustainden. Eine dynamische Eichung ist auch sonst unter Umstanden méglich. ' 

Die Anzeigen der Gerate fiir den atmospharischen Druck, die als Ausgangslinien dienen kénnen, 
sind zumeist rasch zu erhalten, indem man diese z. B. bei Beginn der Aufnahme oder am Ende 
aufzeichnet. Man weiB auch haufig, wo in den Druckverlaufen, z. B. bei pulsierenden Strahltrieb- 
werken, der atmospharische Druck erreicht wird. 


Abb. 5. Unten: Verlauf des statischen Druckes beim Austritt eines Luftstrahles, gemessen am Abb. 6. Verlauf des vorzugsweisen Pitotdrucks 
Rohrrand. Oben: Verlauf des vorzugsweisen Pitotdruckes in einem Luftstrahl beim Rohrausgang. in einem Luftstrahl beim Rohrausgang. 


—s 
At 
Abb. 7. Oben: Verlauf des statischen Druckes am Rohrrand beim Abb. 8. Oben: Verlauf des Pitotdruckes im Strahl eines pulsieren- 
Auspuff eines pulsierenden Strahltriebwerkes, aufgenommen mit den Strahltriebwerkes, aufgenommen mit einem piezoclektrischen 
einem iiblichen Druckelement. Unten: Verlauf des vorzugsweisen MeBgerat mit halbkugelférmigem vorderem Teil. Unten: Verlauf 
Pitotdrucks, aufgenommen im Strahl mit einem piezoelektrischen des statischen Druckes aufgenommen im Strahl eines pulsierenden 
Mefgerat, das vorne eben war, Strahltriebwerkes. 


5. Aufnahme des Verlaufs des Pitotdruckes in einem Luftstrahl und beim Auspuff eines pul- 
sierenden Strahltriebwerkes. Die Untersuchungen wurden so durchgefiihrt, daB auf der einen Seite 
eines Kugelschiebers Uberdruck herrschte, der langsam kleiner wurde (30 atii maximal). Die Druck- 
wellen breiten sich bei der Drehung des Kugelschiebers in einem 4,60 m langen Rohr von 8 cm 
Durchmesser aus und kommen dann ins Freie. Abb. 5 zeigt unten — es ist die hellere Linie — den 
Verlauf des statischen Drucks, wenn er fast am Ende des Rohres an seinem Rande mit Hilfe eines 
ublichen Druckelementes, eines sogenannten Niederdruckgebers, bestimmt wurde. Zur Aufnahme 
wurde ein Schleifenoszillograph verwendet. 

Oben in Abb. 5 ist im wesentlichen der Verlauf des Pitotdruckes beim Rohrausgang aufgezeich- 
net. Hier wurde ein vorne ebenes MeBgerat benutzt. Die Aufnahmeverfahren bei beiden Kurven 
sind ziemlich verschieden, entsprechend der verschiedenen Art des aufgenommenen Druckes. Die 
Ubereinstimmung ist trotzdem gut, soweit dies méglich sein kann. 

Abb. 6 bringt die Aufnahme vorherrschend des Pitotdrucks mit demselben MeBgerat, wobei 
eine Oszillographenschleife héherer Kigenschwingungszahl verwendet wurde. Man erkennt, daB 
kaum Stérungen vorhanden sind. 

Im folgenden werden einige Aufnahmen gebracht, die beim Auspuff eines pulsierenden Strahl- 
triebwerkes von 14 cm Durchmesser gewonnen wurden. Bei diesen Triebwerken wird am Ende 


1 Siehe I § 5 des in FuBnote 7 von Seite 404 genannten Buches. 
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einer Periode, die etwas mehr als "/59 Sek. dauert, Frischluft in das Rohrende eingezogen. Diese 


wird dann zuerst ausgestoBen, wobei fiir kurze Zeit Uberschallgeschwindigkeiten erreicht werden. 
Dies gibt den steilen Anstieg des Pitotdruckverlaufes. Hierauf werden verbrannte leuchtende 
Gase ausgeschoben. 

Abb. 7 zeigt unten ein Oszillogramm, welches mit einem der beschriebenen piezoelektrischen 
MeBgerate zur Aufnahme des vorzugsweisen Pitotdruckes gewonnen wurde, bei welchem der 
vordere Teil eine Ebene bildete. Oben in Abb. 7 ist der statische Druck am Rohrrand fast am Ende 
des Rohres mit einem iiblichen Druckgeber aufgezeichnet. Auch hier sieht man, daB das MeBgerat 
den Vorgingen zu folgen vermag. Die scharfen Spitzen entsprechen nicht dem wirklichen Stré- 
mungsverlauf, sie rithren von der Aufnahmetechnik her; denn die Eigenfrequenz der verwendeten 
Schleife war nicht sehr hoch. 

Abb. 8 zeigt oben dieselbe grundsatzliche Anordnung zur Bestimmung des Pitotdruckes. Doch 


ist der vordere Teil des MeBgerates hier abgerundet mit einer Offnung in der Mitte. Diese Auf- 
nahme wurde mit einer Oszillo- 


graphenschleife hoherer Higen- 
schwingungszahl durchgefiihrt, 
und sie stimmt mit dem wirk- 
lichen Strémungsverlauf iiber- 
ein, der sich aus Flammenautf- 
nahmen ergab. Unten in Abb. 8 
ist der statische Druck mit 
einem MeBgerat aufgenommen, 
das wie sonst bei stationdren 
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Aufnahmen seitliche Schlitze : = 
B = : P INE 
besaB. Erhéhungen des stati- =a = = 
schen Druckes gehen hier nach . ; d 
Abb. 9. Oben: Vorherrschender Pitotdruckverlauf beim Auspuff eines pulsierenden Strahltrieb- 
unten. werkes mit gesteuerten LufteinlaBklappen. Unten: Verlauf des statischen Druckes in der 


Brennkammer. 


In Abb. 9 oben ist der Ver- 
lauf des vorherrschenden Pitot- 
druckes aufgezeichnet am Auspuff eines Triebwerkes mit einem MefBgerat, das vorne eben war. 
Dieses Triebwerk arbeitete mit gesteuerten Lufteinlafklappen. Unten ist der Verlauf des stati- 
schen Druckes in der Brennkam mer zu sehen. 

Wenn auch nur wenige Aufnahmen noch vorliegen, so zeigen sie doch, dafs MeBgerate dieser 


Art gut brauchbar sind. 
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Abb. 10. Temperatureinwirkung bei einem Mefgerat zur Bestimmung des vorzugsweisen Pitotdruckverlaufes im Strahl 
mit ebener Vorderflache. 


6. Temperatureinwirkungen. In den Strahlen der verschiedenen Strahlantriebe kénnen ziemlich 
hohe Temperaturen auftreten. Bei Strahlturbinen sind diese etwa um 600° C, bei Verwendung der 
Nachverbrennung wesentlich héher. Bei pulsierenden Strahltriebwerken betragen die maximalen 
Temperaturen rund 1500° C. Hier strémt jedoch bei jeder Periode zweimal Luft von fast gewohn- 
licher Temperatur an den MeBgeraten voriber. Bei Staustrahltriehwerken und Raketen sind die 
Temperaturen im allgemeinen héher. , wy ; 

Abb. 10 zeigt eine Aufnahme des Pitotdruckverlaufes im wesentlichen mit einem Mefgerat, das 
vorne eben war. Durch das Rohrchen in der Réhrenfeder bis fast vorne in den AbschluBteller hinein 
dringt Kiihlwasser. Man sieht deutlich aus der Aufnahme, die kurz nach Beginn des Laufes eines 
pulsierenden Strahltriebwerkes gewonnen wurde, daB sich die Aufzeichnungen der einzelnen Peri- 
oden kaum andern. Dagegen verschiebt sich die Linic auf dem Oszillogramm, durch welche Atmo- 
spharendruck angezeigt wird. Um das zu verhindern, konnte man mehr Kihlwasser durchflieBen 
lassen, auch die Verstarkereinstellung regeln. Man kénnte auch Kihlwasser durch die auBere 
Geratewand fihren. Doch geniigt im allgemeinen eine Aufnahme iber kurze Zeit. 
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Abb. 11 bringt oben die Aufzeichnung des Pitotdruckes mit einem MeBgerat, das vorne halb- 
kugelférmig war und eine Durchbohrung besaf. Hier sind die Einwirkungen geringer. Unten ist 
die Aufnahme des statischen Druckes mit einem MeBgerat zu sehen, dessen Vorderteil verlangert 
war, wie tiblich, und das drei seitliche Schlitze hatte, durch welche 
der statische Druck wirksam war. Erhéhungen dieses Druckes sind 


i ebenfalls nach unten aufgezeichnet. 


7. Aufzeichnung des verinderlichen statischen Druckes. Zusam- 
= men mit den Aufzeichnungen des Pitotdruckverlaufes wurden schon 
: Aufzeichnungen gebracht, welche den veranderlichen statischen 
= Druck zeigten. Die MeBgerate dazu sehen duferlich ganz ahnlich 

= aus wie ein Hakenrohr fiir stationare Strémungen und sind wirksam 
wie die seither beschriebenen piezoelektrischen Gerate!. Es sind 
Abb. 11. Oben: Temperatureinwirkung Wieder Doppelquarze vorhanden. Senkrecht zur Achse des Gerates 
bei einem Mefgerat zur Bestimmung des sind’ Schlitze oder Anbohrungen angebracht, soweit von der Spitze 


Pitotdruckverlaufes im Strahl mit halb- : ; ‘ : 
kugelférmigem Vorderteil. Unten: Tem- entfernt, daB diese keinen merklichen Einflu8 mehr auf die Auf- 


peratureinwirkung bei einem Mefgerat ° 
zur Bestimmung des statischen Druckes nahme des statischen Druckes haben kann. ; 
Beas Strata seat e ee een en ince ez Abb. 12 zeigt den Aufbau eines MeBgerates. Es wird wieder eine 
nach einer kegelférmigen Spitze. ‘ ig : . 
Réhrenfeder 2 mit AbschluBteller 3 beniitzt, der einen hohlen zylind- 
rischen Doppelquarz gegen ein Widerlager 
7 é I iG 70 driickt. Teil 10 ist eine Spitze. Die Schlitze 
(ZAVLLLLAL LL AEEENSYSY) ii i i A 5 
Kom eee a) eee PS 13 miissen so breit sein, daB rasche Anderun 
y gen des statischen Druckes zum AbschluBb- 
teller 3 zu dringen vermégen. Dazu wurden 
Untersuchungen durchgefiihrt. So haben 
8 8 
sich etwa 3 Schlitze 13 mit Zwischenstiicken 
geringer Ausdehnung und einer Breite von 
Abb. 12. Bebo zur Bestimmung des veranderlichen statischen Drucks. 1mm oder etwas weniger, als brauchbar er- 
I Hohlzylinderdoppelquarz, 2 Réhrenfeder, 3 AbschluBteller, ; 
10 Vorderteil, 13 seitliche Schlitze. wiesen. 
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8. Geriite geringeren duferen Durchmessers ohne Kiihlung. Man kann die Kiihlung weglassen 
und MeBgerate bauen, die als piezoelektrische Teile stabférmige Gebilde besitzen, etwa solche von 
1mm Dicke, 2mm Breite und einigen cm Lange. Man erhalt so MeBgerate geringeren Durch- 
messers. Auch hier gibt es verschiedene Konstruktionsméglichkeiten. Man kann eine diinne 
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Abb. 13. MeBgerat zur Bestimmung der veranderlichen Pitotdrucke oder statischen Drucke unter Verwendung eines Quarzstabes. 
Die Abschlu8Bmembran rechts wird als Vorspanneinrichtung verwendet. 


1 stabférmiger Quarz, 4 Widerlager, 5 Vorspannschraube, 8 elastische Membran, 9 Gehause, 10 Vorderteil. 


elastische Membran als AbschluB und gleichzeitig als Vorspanneinrichtung beniitzen. Dies hat den 
Nachteil, dafs sie Temperaturschwankungen ausgesetzt ist und so ihre Vorspannung andert. Doch 
ist es méglich, solche MeBanordnungen fiir Gase mit geringen Temperaturschwankungen zu ver- 
wenden. Indessen ist es auch vorteilhaft, hier zu Anordnungen iiberzugehen, bei denen die Vor- 
spannung ahnlich erzeugt wird, wie bei den schon beschriebenen MeBgeraten. 


1 F. Staab, Gerat zur Messung des statischen Drucks in einer pulsierenden Strémung, FKFS-Bericht 1944, 
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In Abb. 13 ist 1 ein stabformiger Quarz; 8 ist eine elastische Membran 
Vorderteil 10 des MeBgerates auf das Gehduse 9 an 
verschiebbares Widerlager 4 und eine Schraube 
Leitung zum Verstarker, die mit einer Belegung 
ist mit Hilfe des Leitungsteiles 14 geerdet. 


Die Empfindlichkeit eines solchen Quarzstabes ergibt sich wieder aus der Gréfe des transver- 
salen piezoelektrischen Effektes nach (1). Verwendet man denselben wirksamen Druck, dann erhalt 
man gleiche Elektrizitatsmengen, wenn das Verhaltnis der Ausdehnungen des (rceeribes in y- 
und x-Richtung gleich sind. Dieses Verhaltnis kann ohne Schwiewekeit etwas gréBer genommen 
werden. Doch ist zu beachten, daB® im allgemeinen der wirksame Druck auf die elastische Membran 
und auf den Quarz etwas geringer ist. 

Es ist auch hier méglich, diese Anordnungen zur Bestimmung des veranderlichen statischen 
Druckes im Strahl zu verwenden. Man muB8 nur in dem MeBgerat einen entsprechend ausgebildeten 


Kopf einschrauben. Gerate dieser Art kénnen bis zu einem Durchmesser von einigen wenigen mm 
ausgefiihrt werden. 


welche wieder mit dem 
gepreBt wird. Ein wenig in Achsenxichtung 
5 bilden die Vorspanneinrichtung, und 12 ist die 
des Quarzes I verbunden ist. Die andere Belegung 
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Abb. 14. MeBgerat zur Bestimmung der veranderlichen Pitotdrucke oder statischen Drucke mit einer Vorspanneinrichtung links. 
1 zwei stabférmige Quarze, 2 Vorspannmetallbander, 3 AbschluBstiick, 4 Widerlager, 5 Vorspanneinrichtung. 


Abb. 14 zeigt noch eine Konstruktion, die prinzipiell ahnlich ist wie die seither beschriebenen. 
Diese Gerate kénnen ohne Wasserkiihlung oder auch mit einer solchen ausgefiihrt werden. Es 
sind I zwei in kleinem Abstand parallel gestellte stabférmige Quarze, an deren inneren Seitenflachen 
ein- oder zwei diinne Metallbander 2 anliegen, die zur Vorspannung der Quarze dienen. Diese 
werden auf das Lager 4 gedriickt durch eine Vorspanneinrichtung 5, die verschieden ausgebildet 
sein kann. Die weitere Anordnung ist den seitherigen gleich. Die Quarze I sitzen in einem Ab- 
schluBstiick 3, mit welchem eine ausgegliihte diinne Membran fest verbunden ist. 

Hier ist eine Anordnung ohne Kiihlung eingezeichnet. Es ist jedoch auch méglich, zwischen den 
Metallbandern 2 eine umgebogene Kihlleitung einzuschieben oder anstatt der Seitenbander eine 
elastische Kiihlleitung zu verwenden. Es ist auch méglich, durch den auBeren Gehauseteil Kuhl- 

_wasser zu leiten. Man kann auch etwa drei Quarzstabe beniitzen und diese in der Art eines gleich- 
seitigen dreiseitigen Hohlprismas anordnen. 

Es ist auBerdem méglich, diese verschiedenen Konstruktionen den unterschiedlichen Anwen- 
dungen anzupassen, ebenso entsprechende Verstarker und geeignete Aufzeichnungsapparate zu 
verwenden. 


9. Staurohr. Man kann weiter ein MeBgerat bauen, das gleichzeitig Pitotdruck und statischen 
Druck aufnimmt und als Ergebnis die Differenz, den dynamischen Druck, liefert. Dazu ist es nétig, 
wie im stationaren Fall von Pitotdruck den statischen Druck abzuziehen. 

Abb. 15 bringt die Anordnung. Rechts ist das Gerat, welches den Pitotdruck aufnimmt, weiter 
links das Gerat zur Aufnahme des statischen Druckes mit den seitlichen Schlitzen. Man kann nun 
zwei Quarzanordnungen verschiedener Polung mit gleicher Empfindlichkeit fiir Druckbeanspru- 
chung verwenden und so verfahren, daB man etwa die auBen entstehenden positiven Ladungen 
des Quarzes links zu den negativen des Quarzes rechts hinzufiigt. Anstatt der Quarzstabe kann man 
auch Hohlzylinderquarze anwenden. 
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Da es umstandlich sein kann, zwei Quarzelemente mit einer gewissen Genauigkeit so abzu- 
gleichen, daB ihre Kichkurven entgegengesetzt gleich sind und gleich bleiben, kann man die Diffe- 
renzbildung der entstehenden elektrischen Ladungen auch elektrisch vornehmen. Man kénnte die 
beiden entgegengesetzt gepolten Quarzelemente tiber einen Differentialkondensator an einen Ver- 
starker anschlieBen und ihre verschiedene Empfindlichkeit durch Einstellung des Differential- 


kondensators ausgleichen. Es ist weiter méglich, zwei getrennte Eingangsréhren zu beniitzen und 
einen gemeinsamen Anodenwider- 


stand. Die Empfindlichkeit der bei- 
den Quarze kann hier durch zwei ge- 
trennte Kapazitaten unabhangig von- 
einander abgeglichen werden. AuBer- 
dem ist es méglich, zwei Quarze 
gleicher Polung zu beniitzen und die 
Differenzbildung durch eine Umkehr- 
rohre hervorzurufen, die hinter einer 
der beiden Eingangsréhren Verwen- 
dung findet. Ebenso kénnen Réhren mit entgegengesetzt verlaufenden Kennlinien Anwendung 
finden. Man sieht, auch hier ist es méglich weitgehend den Anforderungen Rechnung zu 
tragen. 


/ Autnahme des 
vA statischen Druckes 


Zh 


——S=S= 4, ; VZ Autrahme des 
SS = fff), Pitotdruckes 


Abb. 15. Piezoelektrisches Staurohr. 


10. Mefgerate zur Bestimmung des verinderlichen Impulses. Bei steil ansteigenden Pulsati- 
onen ist es haufig schwierig, mit Dynamometern etwa den gesamten mittleren durch ein Strahltrieb- 
werk erzeugten Schub im Stand zu bestimmen. Daher liegt der Gedanke nahe, diesen direkt aus 
einer Impulsmessung im Strahl zu gewinnen. In jedem kleinen Zeitelement ist die vom Strahl auf 
den Motor im Stand ausgeiibte Kraft gleich dem Impuls 
der austretenden Gase. Den entstehenden Impuls in 
Achsenrichtung des Strahltriebwerkes kann man _ ver- 
nichten, indem man die Gase senkrecht zu einer Prall- 
platte ablenkt und nach auf en strémen lat. Voraus- 
setzungen fiir eine solche Messung sind, daB sich die Ge- 
schwindigkeiten der Gase von der Diisenmiindung des 
Strahltriebwerkes bis zum Auftreffen auf die Prallplatte 
durch Einwirkungen der umgebenden Luft nicht andern. 
AuBerdem dirfen keine nennenswerten Druckunterschiede 
zwischen Vorder- und Rickseite der Prallplatte entstehen. 


ili Mile Es ist dazu unter Umstanden zu beachten, daB der sta- 
| Sea ees re |e tische Druck in der Diisenmiindung héher als der atmo- 
suk LZ, deel 7 ae : 

4 UL spharische Druck sein kann. 


Es wurden Versuche mit einer glatten Prallplatte durch- 
gefihrt. Doch kann dann eine gréBere Ablenkung des 
Strahls als 90° auftreten. Abhilfe davon wird geschaffen, 
wenn man auf die Prallplatte rund 10 Lagen von Netz- 
geflecht mit einer Maschenweite von einigen mm auflegt. 
Bei Versuchen dieser Art, bei denen aus einem Rohr auf 
eine Prallplatte ein Strahl auftraf, erhielt man bei statio- 
narer Strémung genaue Ubereinstimmung zwischen der be- 
rechneten BewegungsgréBe und dem gemessenen Schub. 


Abb, 16. MeBeinrichtung zur Bestimmung des veran- 


derlichen Impulses im Strahl einer pulsierenden Stré- 
mung. I Prallplatte, 2 Netzbelege, 3 Quarzelemente, 
4 Vorspannhiilse, 7 und 8 Teile der Strémungs- 
leitungseinrichtung. 


Ebenso bekam man geniigend genaue Mittelwerte bei 
Pulsationen nicht zu hoher Frequenz. Der gemessene 
Schub vermehrt sich dazu kaum, wenn man den Abstand 


zwischen Miindung und Prallplatte betrachtlich vergréBert!. 

Bei einem Strahl etwa von 50 cm Durchmesser kommt es darauf an, die Gesamtstrémung so 

wenig als méglich zu stéren. Daher ist es vorteilhaft, einen Strahl geringeren Durchmessers aus 

dem Gesamtstrahl herauszunehmen. Wertvoll ist es dazu, die Kraftaufnahme so durchzufiihren, 

da®B in jedem Augenblick der veranderliche Impuls angezeigt wird. Messungen dieser Art kénnen 

dazu iiber die gesamte Diisenmiindungsflache durchgefiihrt werden, indem man das Gerat in ver- 
schiedenen Abstanden zur Achse aufstellt. 


1 Siehe FuBnote 9 von Seite 404, 
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Zu diesem Zweck wurde in Abb. 16 eine Prallplatte 1 der angedeuteten Art mit Netzbelegen 2 
auf ein Quarzelement 3 gesetzt und fest verbunden!. Das Quarzelement besteht aus einer Vor- 
spannhilse 4, die in titblicher Weise einzelne Quarzplattchen enthalt. Dieser Teil des MeBelementes 
wurde aus der sonst dazu beniitzten Hiille eines vorhandenen Druckmefgeriates herausgenommen, 
welche noch die Kiuhlwasserzufiihrung tragt und das Gewinde, um das Element in einen Motor 
einzuschrauben?. Es war also unnétig das Quarzelement besonders herzustellen. Ferner stellt 5 
das Kabel zum Verstarker dar. Mit Haltebandern 6, die zur Strahlachse parallel verlaufen, wird die 
ganze Anordnung gehalten. 

Vor der Prallplatte ist eine durchbohrte, kegelartige Einrichtung 7 vorhanden, um einen Teil 
der Strémung herauszunehmen. Auf diese Weise wird ein Strahlteil bestimmten Durchmessers aus 
dem Gesamtstrahl des Strahltriebwerkes herausgeschnitten. Der hier gezeichnete scharfe Teil 8 
kann auch unter Umstanden durch einen abgerundeten Vorderteil ersetzt werden. Mit dieser An- 
ordnung 7 wird das Ausstrémen der Gase parallel zur Prallplatte durch die Netze 2 und senkrecht 
zur Achse geférdert. 

Die Vorteile dieser Anordnung sind, daB man auf diese Weise den zeitlich abhangigen Impuls 
des Strahles aufzeichnen kann, auch unter Beriicksichtigung dessen geringer Abhangigkeit vom 
_ Achsenabstand. Man vermag daraus weiter rasch den Mittelschub zu bestimmen. Eine Eichung 
mit Gewichten gibt dazu bei nicht zu hohen Frequenzen brauchbare Werte. Die Stérung der Gas- 
strémung durch das MeBgerat kann abgeschatzt werden. 

Gerate dieser Art wurden mit wenigen cm Durchmesser des Strémungsquerschnittes bei gréBeren 
Strahltriebwerken mit Erfolg verwendet. Bei gewéhnlichen pulsierenden Strahltriebwerken erhielt 
man gute Ubereinstimmung mit Dynamometerangaben. Ebenso waren diese Anordnungen brauch- 
bar bei rascheren Verbrennungen. 


11. Weitere Entwicklungen. Fiir die verschiedenen Anwendungen lassen sich diese MeBgerate 
in gentigender Weise verandern und entsprechend einrichten. Man hat gegenwartig auch Ver- 
starkereinrichtungen, die einen héheren Verstérkungsfaktor besitzen® und kann Verstarker den 
Anforderungen entsprechend bauen. 

Ebenso ist es unter Umstanden méglich, anstatt Quarz andere piezoelektrische Stoffe zu ver- 
wenden mit héherem Piezomodul oder gréBerer piezoelektrischer Druckkonstanten. Man kann so 
die Gerate klein bauen und sie fiir sehr hohe Frequenzen anwendbar machen. So hat das neuerdings 
fiir gewisse Zwecke haufig verwendete Bariumtitanat einen wesentlich héheren Piezomodul als 
Quarz, doch auch eine starke Abhangigkeit der Empfindlichkeit von der Temperatur+. Aus diesen 
Griinden scheint das weniger temperaturempfindliche Keramikelement Zirkonat fiir viele Anwen- 
dungen geeigneter zu sein®. Ebenso kann Seignettesalz gelegentlich Anwendung finden. Lithium- 
sulfat besitzt eine noch héhere piezoelektrische Druckkonstante®. 


12. Zusammenfassung. Piezoelektrische MeBgerate der beschriebenen Art zur Bestimmung des 
veranderlichen Pitotdrucks, des statischen Drucks, des dynamischen Drucks und der veranderlichen 
BewegungseréBe bei pulsierenden Strémungen haben sich gut bewahrt, und es erscheint méglich, 
diese noch weiter den verschiedenen Aufgaben anzupassen. 

Die hier beschriebenen MeBanordnungen wurden 1941 bis 1944 im Forschungsinstitut fiir Kraft- 
fahrwesen und Fahrzeugmotoren an der Technischen Hochschule Stuttgart (KES) entwickelt. Sie 
wurden bei Entwicklungsarbeiten einst auf allen Versuchsplatzen fiir pulsierende Strahltriebwerke 
(z. B. V1-Triebwerk) angewendet. 


Fiir die Erméglichung der Neubeschaftigung mit diesen Fragen wird der Deutschen Forschungs- 
gemeinschaft auch an dieser Stelle sehr gedankt. 


(Eingegangen am 21. Januar 1957) 


Anschrift des Verfassers: Dr.-Ing. F. Staab, Reichenbach a. d. Fils (Wirtt.), Ziegelhof 2 


1 F. Staab, Gerat zur Bestimmung des veranderlichen Impulses im Strabl einer pulsierenden Stromung. 
FKFS-Bericht 1943. 

2 Siehe FuBnote 6 von Seite 404. 

3 Dr. Staiger und Mohilo. Stuttgart-Bad Cannstatt. 

4 Brush Electronics Company, Piezotronic Technical Data (1956). 

5 Atlantic Research Corporation 1956. . 

6 Vgl. das in FuBnote 2 von Seite 406 genannte Buch, Kap. II. 
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Beitrag zur Berechnung raumlich gekriimmter Stibe 
nach der Theorie erster Ordnung 


Von R. Trostel 


1. Allgemeines. Das Ziel jeder statischen Berechnung ist die Dimensionierung eines Trag- 
werkes fiir eine gegebene Belastung. Die Grundlage hierfiir bildet die Kenntnis der Spannungen 
und gegebenenfalls auch der Deformationen, deren Extremalwerte man gewissen zulassigen 
Héchstwerten gegeniiberstellt. Einer Dimensionierung muf daher in jedem Falle die Ermittlung 
des Spannungsverlaufes und méglicherweise auch der Deformationen des Tragwerkes vorangehen. 

Es ist einleuchtend, dai die komplizierte Natur der Tragwerke nur eine naherungsweise Be- 
stimmung dieser GréBen zulaBt dergestalt, daB man den zu untersuchenden Gebilden gewisse 
idealisierende Eigenschaften zuschreibt, die eine erhebliche Vereinfachung des mathematischen 
AufWandes erméglichen. So begniigt man sich beispielsweise hinsichtlich der Deformationen mit 
den Formanderungen einer als Systemlinie bezeichneten gedachten Balkenachse, wahrend man 
hinsichtlich der Spannungen vermége geeigneter Naherungsannahmen tiber die Deformationen der 
einzelnen Stabelemente vollkommen ausreichend durch die Kenntnis ihrer statischen Aquivalenten, 
die man als Schnittlasten bezeichnet, informiert ist. Die nachfolgenden Ausfiihrungen werden 
daher die Ermittlung der Schnittlasten sowie der Deformationen der Systemlinie zum Gegenstand 
haben. 

Bei der Bearbeitung statischer Probleme unterscheidet man zwei grundsatzlich verschiedene 
Betrachtungsweisen. Die einfachere setzt die am System auftretenden Deformationen hinsichtlich 
der Lage der Lasten als vernachlassigbar klein gegeniiber den Systemabmessungen voraus, so daB 
sdmtliche Schnittlasten am unverformten System berechnet werden kénnen. Dieser Theorie erster 
Ordnung steht eine genaue Theorie gegeniiber, bei der die Schnittlasten von vornherein fiir das 
verformte System bestimmt werden. Naturgema8 ist diese Problematik wesentlich schwieriger, 
und eine Liésung gelingt nur fir spezielle Falle mit ertraglichem Rechenaufwand. Wir wollen nach- 
folgend nur die Theorie erster Ordnung zum Gegenstand der Untersuchung machen. 

Die maBgebenden Gleichungen fiir eine statische Berechnung sind die aus den allgemeinen 
Bewegungsgleichungen 


ER, =mi,, sm=f |e x tdm 
fiir den Fall beschleunigungsloser Vorgange sich ergebenden Gleichgewichtsbedingungen 


a On, 2 Nie 0 (la, b) 
sowie spezielle Formanderungsrelationen. 
Die mathematische Darstellung der Systemlinie ist mittels eines, von einem skalaren Parameter t 
abhangigen Ortsvektors 


T= T(t) = {a(t)s ¥(0)5 2()} (2) 
méglich. In jedem Punkte einer solchen Raumkurve sind drei zueinander orthogonale Richtungen 
durch die Einheitsvektorent 

oe t(t) _—t(t) dt dt dt/ds 
£ nal 


[EO] [@l dt ds” = Kn Mable b= ¢.= tx il, ) cape) 


definiert. Von diesen interessieren im Hinblick auf die Berechnung der Deformationen noch die 
Ableitungen nach der Bogenlange s, die vermége der Frenetschen Formeln wiederum durch die 
Einheitsvektoren selbst ausdriickbar sind: 


dee = dey 


de 
ds eae ds =—Kes +x ey, —— 


7 ee (4, a, b, c) 


Hierbei bedeuten K = 1/R die Kriimmung und x = da,/ds den Verdrehungswinkel je Langeneinheit 


(Windung), um den sich das durch t, 1 und b gebildete Dreibein beim Fortschreiten auf der Raum- 
kurve um diese dreht. 


* Der Punkt deutet eine Ableitung nach dem Parameter t an, und s ist die Bogenlange der Raumkurve. 
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Als einschrankende Voraussetzungen gelten fiir die nachfolgende Rechnung drei Annahmen: 


cA Systemlinie sei die Schwerlinie und gleichzeitig die Verbindungslinie der Schubmittel- 
punkte. 


2. Die Belastungen sollen unmittelbar in der Systemlinie angreifen. 


ay = Richtungen der Haupttragheitsachsen sollen in jedem Punkte der Systemlinie mit den 

ort vor iegenden durch ey, baw. ee definierten natirlichen Richtungen der Raumkurve identisch 
sein, was immer dann erfiillt ist, wenn eine der Haupttragheitsachsen durch den, dem entsprechen- 
den Punkte der Systemlinie zugehérigen Kriimmungsmittelpunkt verlauft. 


2. Berechnung der Schnittlasten nach der Theorie erster Ordnung. Schnittlasten sind Uber- 
tragungslasten zwischen den einzelnen Teilen einer Konstruktion. Sie werden aus den Gleich- 
gewichtsbedingungen fiir ein abgeschnitten gedachtes Teilstiick eines Systems berechnet. Ent- 
sprechend der beiden Gleichungen (la, b) unterscheidet man (Abb. 1) 

1. Schnittkrafte O(s) bzw. Q(t) und 

2. Schnittmomente Ji(s) bzw. NU(t). 


| (ojo 


q (0), m (0) 


eg (ft) 


< M/s) 
eg (b) igi we -O(s) 
ss <IM(S) a 
ae De) : 


-O 
Ee 0 


-M, 


Abb. 1. Abb. 2: 


Sie treten an zwei sich gegeniiberliegenden Schnittufern nach dem Reaktionsprinzip stets paarweise 
und mit jeweils umgekehrtem Richtungssinn auf. Indem man diese Schnittlastvektoren skalar mit 
den Einheitsvektoren ez, ¢, und e; multipliziert, ergeben sich die Langs- und Querkrafte bzw. 
Torsions- und Biegemomente. 

Das Krafte- und Momenten-Gleichgewicht des in Abb. 2 dargestellten Systemteils, der durch 
Krafte und Momente singular (x, Ix) bzw. kontinuierlich (q(s), m(s)) belastet ist, fordert 


— Ly + & Be + fale) do + Os) = 0, 


also 
Q6) =O), — 3S Be — fale) do, (6) 
— My — (to — t(s)) X Do +3 Mee + J m(a) do +3 (tx —r(s)) X Be 
+ i; (t(o) —1(s)) x g(a) do + WM(s) = 0, 
bzw. 


Ms) = My + (t= Hs)) Do — SMe — J mo) da — ¥ (te —160)) x Be 


— f (t(c) —t(s)) x g(a) do. (6) 


o=0 
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Die Gleichungen (5) und (6) erméglichen die Schnittlastermittlung, sofern die Anfangswerte Qo 
und St, vorgegeben sind. Bei Systemen mit freien Enden wird man daher den anfanglichen Orts- 
vektor t, dem freien Endpunkt zuordnen, da in diesem Falle D, = Yi, = 0 gesetzt werden kénnen. 
Ansonsten sind Q, und S){j meist mit Auflagerlasten identisch, die aus einer Vorberechnung mittels 
der Gleichgewichtsbedingungen fiir das Gesamtsystem hervorgehen, sofern dieses statisch bestimmt 
gestiitzt ist. Im Falle statisch unbestimmter Stiitzung miissen zu den Gleichgewichtsbedingungen 
noch gewisse geometrische Bedingungen der Kontinuitat hinzugenommen werden, wovon am Ende 
der Ausfithrungen noch gesprochen wird. 

Durch Differentiation von (5) und (6) erhalt man 

dy dm dt 


=A), =a A X Als) —M) =— eg x Dls) — ms) (7a, b) 


oder, wenn man diese Differentialbeziehungen hinsichtlich des Parameters t entwickelt, 
Oi) =—li]q@, MO) =— Ie] [ee x O} +m]. (8a, b) 


Bei einer Komponentendarstellung von (7a, b) ist zu beachten, daB bei der Differentiation der 
Schnittlastvektoren, sofern ihre Komponentendarstellung auf die Basisvektoren ¢;, ¢,, und ¢; Bezug 
nimmt, diese mitzudifferenticren sind. Man hat also z. B. wegen 


Y= {023 Q,3 Oc} = Qs es + Q, Cn + Qe ee 


nach Differentiation 


dQ, dOe |, Wr 
ds § ds ' ” ds 


woraus schlieBlich unter Beachtung von (4a, b, c) 


dQ, d d d 


und damit aus (7a) 


dQ¢ deg ! dey, deg 
r ee ds ! 27 Qn as ! Ora 


dQ: dQ dQ¢ 
Tee i es ee KR Oph Or ds qa ty am wt 
hervorgeht. Analog hierzu lautet die Komponentendarstellung von (7b) 
dM: dM, dM; 
= eM ne, Se ae 4, + *M, = —m,—Q,- 


3. Berechnung des Deformationszustandes. Wie bereits betont, begniigt man sich bei den Ver- 
formungen naherungsweise mit der Darstellung des Verschiebungs- und Verdrehungszustandes der 
Systemlinie. Man geht hierbei, ebenso wie im allgemeinen in der Elastizitatstheorie, von den Ver- 
formungen der Elemente, im vorliegenden Falle also der Stabelemente aus. Fiir diese werden fol- 
gende Annahmen gemacht: 


1. Die Elementendeformationen infolge der Schnittmomenten-Belastungen 4uBern sich in 
reinen gegenseitigen Querschnittsverdrehungen ohne Langenanderungen der Systemlinie, und die 
ebenbleibenden Querschnitte! stehen nach der Deformation weiterhin senkrecht zur verformten 
Systemlinie (Abb. 3). 

2. Die Elementendeformationen infolge der Schnittkraft-Belastungen auBern sich in gegen- 
seitigen Parallelverschiebungen benachbarter ebenbleibender! Querschnittsflachen (Abb. 4). 

In der Festigkeitslehre werden die je Langeneinheit auftretenden gegenseitigen reinen Quer- 
schnittsverdrehungen y™)aus der Momentenbelastung und die Parallelverschiebungen ¢@) aus der 
Schnittkraftbelastung durch die Relationen 

\ 
(ez ° ° 


= I 

70) = fy; yl; 00} = Cys Ms er ea | = 2 Oa (9) 
| 1 
(Oe oars, 


1 Vom Einfluf jeglicher Querschnittsverwolbungen bei der Torsion und bei den Querkraftbelastungen 
wird hier also abgesehen. 
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und 


ra(O) : ; _ ; : 
é {e); Ws VO} = {0z; Ons Qi} |)-0 %, GF 0 =D Io (10) 
1 
0 0 cer 


mit den Schnittlasten in Beziehung gesetzt. Hierbei bedeuten J, bzw. J, die Tragheitsmomente 
hinsichtlich der Achsen 4 baw. ¢, J, das Torsionsflachenmoment, dak im wesentlichen durch das 
Volumenintegral des Hiigels der Torsionsfunktion definiert ist und «, bzw. X- zwei die Schub- 
spannungsverteilung in der Querschnittsflache F beriicksichtigende F piegre! 

Aus diesen Schnittlast-Verformungsrelationen sind nunmehr die Differentialgleichungen fiir die 
DeformationsgréBen des Systems zu entwickeln. Der Verformungszustand ist eindeutig beschrieben, 
wenn man die Verschiebungen tv bzw. die Verdrehungen @ der Punkte der Systemlinie, sowie die 
Verdrehungen 6 der Flachenquerschnitte berechnet hat. Wir suchen daher nachfolgend, die ge- 
nannten GréBen mit den jeweils dreikomponentigen Dehnungsvektoren y™) und ¢@) in Verbindung 
zu bringen,. 


Abb. 3. Abb. 4. Abb, 5. 


Vereinfachend wird zu diesem Zwecke der gesamte Deformationszustand in einen Schnitt- 
momenten-Deformationszustand (wv), G™)) und in einen Schnittkraft-Deformationszustand 
(w®, p@, B@) zerlegt. Die jeweils resultierenden Relationen lassen sich superponieren, womit 
dann die endgiiltigen Deformationsgleichungen gefunden sind. 

a) Die Schnittmomenten-Deformationen. Da die Flachenquerschnitte hierbei nach der 
Deformation weiterhin senkrecht auf der deformierten Systemlinie stehen sollen (Abb. 3), sind die 
Querschnittsverdrehungen f\™) mit den Verdrehungen @™) der Systemlinienpunkte identisch. Man 
hat demnach fiir die Verdrehungsanderungen der Systemlinienpunkte die Relation 

a(M M 
Fe = Ee =H) = MN Oy (11) 

b) Die Schnittkraft-Deformationen. Im Gegensatz zu a) sind hier die Komponenten des 
Dehnungsvektors € @ nicht allein durch die VerformungsgréBen der Systemlinie (Ww, gp) 
angebbar, da eine Schubverzerrung der Elemente sowohl eine Neigung g@ = {p@; g); gy} 
= {69; g; gy) }der Systemlinie als auch eine (wegen der Parallelitat der Querschnitte bedingte 


konstante) Verdrehung pe der Querschnitte aus ihrer unverdrehten Lage heraus zur Folge haben 
kann. GemaB Abb. 4 hat man also 


dy(Q) 
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oder auch 


S 7 div(@) - 
GO — BP) x ey + (MO eg) eg = 2 = 0.8, (12) 
und nach vektorischer Multiplikation mit ¢; 
FO — BY) — (oO — BE) ep = — BM = eg x (O99). (12a) 
Durch Differentiation von (12a) gewinnt man 
Fag 
Sr re (OQ Jo)] » (12b) 


so daB nach Zusammenfassung von (11) und (12b) 


Fre OEE d 
F— + = F = MO + 5 [ee x © F0)] (13) 


entsteht. Um hieraus auf die Differentialgleichungen fiir die Verschiebungen zu kommen, betrach- 
: dw 
ten wir noch einmal ein Element (Abb. 5), dessen relative Verschiebung dy = oe ds des Punktes B 


gegeniiber A man sich hervorgegangen denken kann einerseits durch eine reine Streckung 
e®) ds = de; des Elementes und andererseits durch eine Drehung um A um den Winkel 9. 


Man hat demnach mit. 


d3' = (1 + &:)d38 ~1-d3 =dse, 
die Beziehung 
= d x 
dw = (dw e,)ee + y X 43 = (2 as es) eg + ds are 
also 


d\ d\ = Q = 
% = (gees) + OX = (Beles +B X & (14a) 


bzw. nach vektorischer Multiplikation mit Cs 


dw - 
CE yg PSE (14b) 
und damit nach Differentiation 
ded d dw 
To ge ts ea (e x @). (4c) 
Einsetzen von (14c) in (13) liefert damit schlieBlich 
d d div d 
gs Pe &s) Pg, (e x =) = MO + z, [ee X (Q99)], (15) 
also in Komponentenschreibweise unter Beachtung von (4a, b, c) 
dye dwe _ Ms Qc 
ae (a x, EJ, 0s OR (15a) 
SN, pe fae _ M, Qn d/ @ 
ane (a xy) “(Te a Kw; — xp] = Fee oe ie, Galaakads Exaok Oe) 
d (dw dwe wugeay Qe d /_ Qn 
ds & + K te — 7 ue) — 0 ( ie 7w,) sR fee Li yim ane (ser): (15¢) 


Zu diesen drei Gleichungen tritt nun noch die bisher nicht ausgenutzte ¢-Komponente der Glei- 
chung (14a) 
dio dwe preanUe 
abo ok ad ee 


hinzu, so daB mit (15a—d) vier Differentialgleichungen fiir die Verschiebungen 1) = {we; w,; we} 
und den Drillwinkel pm; zur Verfiigung stehen. Nach Lésung ergibt sich fiir die beiden anderen 
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Drehwinkelkomponenten gema8 Gleichung (14b) mit (15a) und (15b) 


_ __ (dwg _ 1 (dpe Me Qe 
Pa (z + %0,) rate ag acG F’ 
loa—b) 
dw 1 d/ldgp:\ . K l1d/1M. 1 M. Q ( 
ee Rap ye, = ae es af a n 
eee gee eae 7 as (R i) Warp ra (K ES) EI, ) a, GE” 


wahrend die resultierenden Querschnittsverdrehungen f sich unter Beachtung von (12a) zu 


B = BOO + BR = 9 + BY = GM) + HO) — ey x (O9) =F—ee x OD), — (17) 
also in Komponentenschreibweise zu 
: ws oat Q Q 
Be ee ce = Pn a nF 2 Be a = Er (18a—c) 


77) 


ergeben. 
Durch fortlaufende Elimination laBt sich das Gleichungssystem (15a—d) auf zwei Differential- 
gleichungen dritter Ordnung fiir y; baw. w; reduzieren. Zunachst erhalt man aus (15c) nach 


Elimination der Verschiebungskomponenten mit Hilfe von (15a) und (15b) eine Differential- 
gleichung dritter Ordnung fiir den Drillwinkel ¢,: 


d 1d 1 dpg\| , d /K : x dpe _ x Mz dads AM: aes Mz 19) 
ds |x ds\K ds ' ds reg ' K ds KG Ji mG ds (x &3,)| ! ds (7 Ey) Dei bts Fe ( 
Nach Lésung von (19) kann nunmehr aus (15a) und (15b) unter Eleranziehung von (15d) folgendes 


Gleichungssystem fiir w, und w; gewonnen werden: 
dwe “udwe 1 Me , *% Qe , Qe 1 dye 
ds K ds KG Ii KE Fa” apG.h aK ods? 
Ky, 4 af Ldwe) 1 d/l Qe) 1 pt d/l doe 
bere 8 oe de\K ds % ds (k EF) «|x ds\K ds 


Becks 1 d/l. M; 1 M, Qn 
o nz PS x ds (Kea) ~ % EJ» t CF: (20) 


(20a) 


Damit folgt schlieBlich durch Elimination von w, gema (20b) aus (20a) unter Beriicksichtigung 
von (19) eine Differentialgleichung dritter Ordnung fiir w;: 


d{1 d/l dw¢ d /K , % dwe dilgland sla O: 
ds |\% ds\K ds}| ' ds\x “§)" K ds ~ ds vas (KEF 
1 Mz _dge\ | K Mt , K d/ Qn Qe Qe 
| K|2 (es, le * EJ: | va (,0r) + “art Kaen 


1 d/l dg\ , K ede IMME oo ol ee On aed: fl (21) 
wz ds\K ds} ' % %6 walKeH| % EJ, | %GF\ds\x}° 


Nach Lésung des Gleichungssystems (15a—d) stehen also insgesamt 6 Integrationskonstanten 
zur Befriedigung des Randwertproblems zur Verfiigung. Lassen sich mit diesen Konstanten die 
aus der Stablagerung hervorgehenden geometrischen Bedingungen nicht erfiillen, so liegt statisch 
unbestimmte Stiitzung vor. In diesen Fallen bleiben bei der Schnittlastenermittlung nach (5) 
bzw. (6) und damit auch in den Lésungen von (15a—d) so viele Fessellastkomponenten (statisch 
unbestimmte GréBen) zunachst willkiirlich, wie notwendig sind, um zusammen mit den 6 Inte- 
grationskonstanten das Randwertproblem zu erfiillen. 

Vor der Untersuchung eines Beispieles seien noch die Gleichungen fiir Stabe mit ebener Anfangs- 
kriimmung angegeben. Mit x = 0 gehen aus (15a—d) die Beziehungen 


dye a Me Qc 


ds | ds GJt aeG F’ 
d?we M, df On } 
Reg ny, § ds lec P| ; 


(22a—d) 
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sowie aus (16a, b) die Beziehungen 


dwe dw an 
0, == 53 (= 7 + Kus (23a—b) 


hervor. Auch hier ist eine Reduktion des Systems (22a—d) auf zwei Differentialgleichungen dritter 
Ordnung méglich. Setzt man die aus (22d) hervorgehende Relation 


Poe dwg Qe 24 
be Bama (a EF (24a) 
in (22c) ein, so folgt eine Differentialgleichung dritter Ordnung fiir we: 
Pup dened Me. paces 0, 
as (K Z| gs WA #2) = Fey, ds lac VK Eh) ee on ce) 


wahrend man nach Elimination des Drillwinkels gy; gemaB (22b) 


bee eee, chi Qe 24, 
V5 = Reise eR, By, fe Reds aOR) ee 


aus (22a) die Gleichung 
d[l@ue] , p-dw;_ Ms d/l M, 
ds IK ds ds GJ: ds (x Ed) 


Q@ , diid/ & ; 
aeGE ' ds Ik ds a a] a) 


erhalt. 


4. Beispiel. Eine zylindrische Schraubenfeder mit n Win- 
dungen, dem Windungshalbmesser r und dem Drahthalbmesser a 
besitze an ihren Enden je einen starren, bis in die Zylinder- 


achse reichenden horizontalen Arm OA bzw. O'A' (Abb. 6). Es 
soll der Deformationszustand der Feder infolge zweier, an den 
Armen angreifender Krafte P bestimmt werden unter der Vor- 
aussetzung, daf} der Kraftangriffspunkt 0 als unverschieblich 
angesehen wird. 

Die Gleichung einer Schraubenlinie vom Zylinderhalbmesser r 
und der Ganghéhe h lJautet fiir das in Abb. 6 eingetragene Ko- 


ordinatensystem mit 


s(t) =Verr ae l= t/t ae (e-) = 4 eyt ros, 


296 
i yh 
POS oon? 
Ss 
alsoute—t(s)—— folgendermaBen: 
(s) Ti Tee 
C= 7 ,C0s t; sink; oT} 
Ss S Ss 
= Peau ey ten 3. ‘sin ; a ———— 
| rJl+e : ryl4+ a ; rjl+ 
Die Basisvektoren folgen nach (3a, b, c) zu 
dt 1 
nls Speaker a ; 
er ja sint; cost; a}, 
__ deg/ds ; 
& EEREET {cos t; sint; 0}, 
er = es XC : {x sin t 1 
== les SS ee OS LI eX COSINS 
¢ 5 n yi ag ee ? a 


und fiir Kriimmung und Windung erhalt man aus (4a, c) 


1 1 y 


K=—= ay en ae iret 
R r (1 + a) Ky R,’ 4 = (ee aan meer 
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Die Schnittlasten ergeben sich nach (5) und (6) Mamit tg 03, = 0, 0, = = PK0; 07 1), 
q(s) = m(s) = 0, By = 0, MN, =02 


As) = Dy = B = P {0;0; 1} 


und 
Ms) = — r(s) xX Oy = —1(s) oe Pr {—sint; cost; 0}, 


so daB man fiir die Komponenten 


Q:=De-= P. = 5 Oe. 0; Cen Poe 


1 2 2” 
V 7s Vi+e (25) 
=e Pro; M =e = 0; Me=Me= = 
| Sas pe aS etal ee eee 
x : : ae dt d 1 d ; c a at 
erhalt. Damit ergeben sich aus (19) baw. (21) mit ae de ae tre @ sowie mit J, = ae 
x at F x 3 E : 
ee re =A", O =H%=7> © = 37 Fw Me Gleichungen 
dye dye 4Pr?v @ 
a hb Baa Te 2°) 
; . () =e a = 
bzw. wenn man meee zur Abkiirzung EF = Ewai jie = Ep Setzt, 
Gwe , dwg __ 4 Pr? 21+» /a\2 dg: , a 
ee eg ae lpea(i +27 eer Gtem tr ingens eo) 
Zunachst ist Gleichung (26) zu betrachten. Ihre allgemeine Lisung lautet 
4 Pr . 
G: = C,4 ret Taga | C,cost + C,sint, (28) 


so daB (27) in 
dw: ei 
so de a 


(os 


1 2 2p / 
ala | yo 2 ar (C, sin t — C; cos t) 


mit der allgemeinen Lésung 
4 Pr 2 2 1+» /a\2 CE IPE ; 
= Feat! ad ea ae roel) |! |- =! | BCE SOe Lien Ca ety), 
+ C, + C; cost + C, sint (29) 
iibergeht. Die restlichen Verformungsgréfen ergeben sich dann nach (15d), (20b), (16a, b) und 
(18a—c) zu 


Wide ie 2 21+» /a\2 rE Se 
w, == yl 4 oa I! Tare L rata ria ar (C,tcost + C,t sin t) 


+ (C, + ar C,) cos t — (C; + ar C,) sin ; ; (30a) 


4Pr 1—« 21+» /a\2 Cates C, 1+ 2? 
eos. Poel agree aac Mirage G i a C 


+ o2r (Cyt sin t— Cyt cost) — [x C, + r(1 +a?) C,] cos t-— [a Cy + r (1 +07) Cj] sint, (30b) 


p= — Fo ral Pao sire (£)") + Gsin e— Go0s1], (30c) 
Le z 4 ye o (C, cost + C, sin t), (30d) 
Br = = G4 yea | C, cost + C,sint, (30e) 
r= %+ ap rea afen( — s) + C sin t— C; cos (30f) 
Bas ore by es pata (Cy c08t + Cy sin). (30g) 
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Die Konstanten C,, (k = 1... 6) werden aus den Randbedingungen bestimmt. Sehen wir den 


Kraftangriffspunkt 0 als unverschieblich an, so ist fiir den Arm OA nur eine Drehung um O miglich. 
Der Verbund zwischen Arm und Schraubenfeder bedingt, daB der Drehwinkel des Armes mit dem 


Querschnittdrehwinkel B (t =0) des Schraubenfederpunktes A identisch sein muB. Man hat also 
fiir die Verschiebung des Federpunktes A 
wy, = Ww (t = 0) =P (t=0) x (OA) =—B (¢=0) x re, (t= 0) baw. [w +rB X egh-o=9 (31) 
oder in Komponentendarstellung 

w; (t = 0) —rB,(t=0) =0; w,(t=0)=0; w(t =0) + rf = 0) =0. (3la—c) 


Als weitere Randbedingungen betrachten wir die Aussage, daB die Verschiebung des oberen 
Kraftangriffspunktes O' nur eine lotrechte Komponente tv,/ = w,/¢, besitzen soll. In diesem 


Falle folgt mit dem Querschnitts-Neigungswinkel fa (t =2 an) (der mit dem Neigungswinkel des 


Armes 0’ A’ identisch ist) fiir die Verschiebungen des oberen Schraubenfederpunktes A’ 


wa = w(t =2an) = 1, —r f(t =2 27) x @, (t =2an) baw. [w a a Cnlia2an = Wey 


(32) 
also in Komponentendarstellung 
IN Pi Jer ESCO, OP ae w, (t =2an)=0, 
; (32a—c) 
we (t = 2707) SN ae A icra eres : 
a 


Einsetzen von (29), (30a), (30b), (30e), (30f) und (30g) in (3la—c) baw. (32a—c) liefert ein lineares 
Gleichungssystem fiir die C, (k = 1,...6) sowie fiir w,’. Hieraus folgen schlieflich 


r 


C= 0 >Co— 0, Cs Oi 
ee eee 
wg = Barns, +2282 Falls at DLEB EY 
erat tat teas dey), a 
wobei C, zunachst unbestimmt bleibt. Die Deformationsgré&en sind also endgiiltig 
ah A Se a - | (34a) 
is pealite|l +ryat areal) |: - (34b) 
ele erate (34) 
mec eee (34a) 
pS GA 
0% Ena l+q@” (34e) 
wear Ene |! bs | 2 iFal)|¢ sin t) 4 Fear ptt t sing), (34f) 
ae en aba 
ory eo ee eel ; marci (1 — cos #) , (34g) 
BSG ees (Ue Oca eee 


2» 21l+y7/a 2 
1 F a 1 
>| Depress rpa(s) est toe pear) @— sing], (34h) 
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Zur Deutung der in der Lésung verbliebenen Konstanten C, betrachten wir ein Element der 
Schraubenlinie (Abb. 7) mit den Verschiebungen w; = r C,/x und We = —rC,. Man erkennt, daB 
der hieraus resultierende Verschiebungszustand, da sich fiir den Steigungswinkel w der Schrauben- 
linie tg @ = h/27r =o ergibt, einer waagerechten Verschiebung des betrachteten Schrauben- 
linienpunktes entspricht. Durch C, wird also den Verschiebungen samtlicher Schraubenlinien- 


punkte eine konstante, horizontale, in die Tangentialebene an den der Schraubenlinie einbeschrie- 


benen Zylinder fallende Verschiebung 


Case ule 
wp = — V1 +03 ie 


iiberlagert, die nur durch eine verformungslose Drehung 
der gesamten Feder um den Winkel 


wT 


ae 
ayia rc 
Le ae Re yi +o, Bea Iya 


also durch die Drehungen Abb. 7. 


C. 
Pe = Pz hs lg = C, und Y alia EG Whe a =f 


[vgl. auch (34a) und (34c)] verwirklicht werden kann. Durch C, wird also eine beliebige verfor- 
mungslose Gesamtdrehung des Systems um die z-Achse beschrieben, so da wir C, = 0 setzen 
kénnen, wenn wir eine solche Drehung ausschlieBen wollen. 

Mit ¢ w,/ = P folgt aus (33) fiir die Federkonstante einer Schraubenfeder 


P E a‘ 1 
(vss Wz/ = or Gaaer 14 v a 21+» ey ] 6E (2) . (34) 
at YT seve Sale eeeT 4 1+?\r 


Sie geht fiir den Fall sehr flacher Windungen (« = h/2 ar < 1) und fir im Vergleich zum Win- 
dungshalbmesser r kleine Drahthalbmesser a in die bekannte Form! 


Ea G at 


o> 8rn(Q+v) 4rn 


iiber. 


(Eingegangen am 13. Februar 1957) 


Anschrift des Verfassers: Priv.-Doz. Dr.-Ing. Rudolf Trostel, Lehrstuhl fiir Mechanik (Prof. Dr.-Ing. Szabo) 
der Technischen Universitat, Berlin-Charlottenburg 2, HardenbergstraBe 34 


1S, z. B. I. Szabé: Einfiihrung in die Technische Mechanik, 2. Aufl. S. 137, Berlin/Gottingen/Heidelberg 
1956. 
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Free and Forced Vibrations of Straight Elastic Bars According to the 
‘(Method of Internal Constraints” 


By E. Volterra and E. C. Zachmanoglou 


1. Introduction. Professor S. P. Timoshenko in two well-known papers! has discussed the 
effects of transverse shear and of rotatory inertia on the natural frequency of lateral vibrations 
of a simply supported bar. He has shown that the correction for shear in the frequency equation 
is, in this particular case, about four times that due the rotatory inertia of the bar. The equation of 
motion first derived by Timoshenko has been used by E. Goens? to discuss lateral vibrations of a 
free-free beam, by R. M. Davies® to discuss lateral vibrations of a loaded fixed-free bar and by 
E. T. Kruszewski4 to discuss lateral vibrations of a cantilever and of a free-free beam. R. W. Traill- 
Nash and A. R. Collar® have applied Timoshenko’s analysis to derive frequency equations and 
modal functions for flexural vibrations of uniform beams in the more general cases of end con- 
ditions which can be derived by the various combinations of the three usual types of end supports. 

In recent papers® the equations of motion for curved bars and curved and twisted bars 
have been derived taking into account the effects of shear and of rotatory inertia by applying the 
“Method of Internal Constraints“. It has been shown how these equations can be solved in terms 
of the eigenvalues and of the eigenfunctions of the straight bars cases. It is the object of the present 
paper to discuss the problems of free and forced vibrations of straight bars of uniform cross-section 
taking into account the effects of shear and of rotatory inertia by applying the “Method of Internal 
Constraints“‘. 

The equations of motion for free vibrations derived in previous papers® for straight bars can 
be written under the following forms: 


SS 6) Bay, or Oy 5 Ax 
(l-20)(1—20) “ext Ci Oe ore 
a ©) 
E aes et) — 0 2 
ax? | Ox a” 
(—o)E  .0us Oe 9 Olly 
(Geo) (22 ee re 
ES i LE oe o 
Ox? ax & ae? 
aA ov OPA 
2 None ee ee 7 se y 
C raga — A+ 26)4,—1 FT A w= ert ae? 
2 Or Uz OV 4 Cruz 
GC ya — A+ 26) w,—AZ Ah, = 075 or ’ (3) 
OVx, Ohy Oz OV 
(20) ra Spee oe a ors | 
OAs a Ons 
Grizg—Gi,_ Cu, =0t: =. 
4) 
2 Ply Oty 
Cry a Gu, —GA, = ly ae 
with the following boundary conditions atx = 0, x = L: 
Ohx Ovy 
EB Wei. E 3 A. bv, =0, (5) 
* S. P. Timoshenko, Philos. Mag. 41 (1921) p. 744 and 43 (1923) p. 125. 
® E. Goens, Ann. Phys. 11 (1931) p. 649. 
° R. M. Davies, Philos. Mag. 23 (1937) p. 653; id id id (id) p. 1129. 
3 diy, TE, Kruszewski, N. A. C. A.; T. N. 1909 (1949). 
: a ve Traill-Nash and A, R. Collar, Quarterly Journal of Mech. and Appl. Math. 6 (1953) p. 186. 
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Ox 0 z 
on Ox OV, 
fz] —0, [2] 5. =0, (+26) 32424, +Ap,|d0, =0, (7) 
Ou Oh: 
A Ofly =a | OA, — 0h ‘e) 


The above equations have been derived by assuming that during motion the components u, v, w 
of the elastic displacements are expressed under the following form: 


u(x, Y> z3 t) = V(x, t) =p vAAx, t) oe z [xl X, t) Z 
v(x, ¥, 25 t) = v,(x, t) + yA,(x, t) + 2 p(x, t), (9) 
w(x, ¥> &; t) a ¥,(x, t) an Of A(x, t) as od [U(X t) 


Vax Vy> Ves Ags Ayr Aer Uxs fly, and yt, being nine unknown functions of the variable x and of the 
time t, and by the application of Hamilton’s Principle. The bar is referred to a system of ortho- 
gonal cartesian coordinates 0, x, y, 2, where x coincides with the axis of the bar, while y and z 
with the principal axes of inertia of the cross-section of the bar, A is Lame’s constant and G is the 
modulus of rigidity of the material of the bar related to Young’s modulus E and to Poisson’s 
ratio o by the relationships: 
Ne ent (Giga. 
(+0) @—20) 2(1 0) 


While r,, r, are the radii of inertia of the cross-section of the bar with respect to the z and y axes 
and @ is the density of the material of the bar. 


In the case of two-dimensional elasticity the equations which determine the four unknown 
functions /,, /,, , and ¥, are 


eA, E er av 
2 Jee _y2 ad ios BAAN, 
CT a ~ A —o)" tx Gi,—C ae a 
Ovy —  Ovy Ode 
i eye 
Oye, E Ove Ey @dhy 
0 a ~ OC —o0%) tx? © (—o) ox’ oo 
poly eae ee Pe Eo a0, 
Oe Ox (1 —o%) “9 (1—o®) Ox 


with the following boundary conditions at x —0, x = L: 


Ax ov 

Hae ae be PAa|Oes = 0 i 
Ox hy ies is 
ae + o%, iv, =0, |, =. (13) 


2. Lateral Eigenvibrations. In order to study lateral eigenvibrations of the bar, the method 
of separation of variables will be used. Let in equations (1), (2), (5) and (6) 


2%, 1) = M(x) AQ > Mal® 1) = M(x) Ai) | 
r(x) = Vx), v(x, 1) = Vix) i) - 


Substituting equations (14) into equations (1) or (2) after changing r, or ry with r, one obtains 


(14) 


(l—o)E PX, (x) «) — Ga) 
fee ee d& _ 1PAO_ __ op 
or? X,(x) OV oe " (15) 
CRIN) dX,(x) 
gan Cc LAC ee 


0 Yi) OM ue 
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and since the constant «, is a real quantity!, lateral eigenvibrations are defined by the following 
equations obtained from equations (15) and the corresponding boundary conditions: 


oe 0 


dt? 
l1—o)E 9 aX. dy. 
al ' 6) we 2 ae a — G X,(x) —G oe + wi or X,(x) =0, (16) 
CPI 6 rote nero | 
Px a0 
‘ atx — 04 — Ee (17) 


The eigenfunctions X, Y, satisfy the following orthogonality property* 
L 
IMOGe ORs ap i Xin Xm) dx =0 for Min - Wim: (18) 
0 

By eliminating X, between the two last equations (16) the following equation in Y, is obtained: 


GING @Y, w@of[E , l+o0)—20) 
ae og “ae E E v (1 — 0) 


(1+ 0)(1—20) | Be =: (19) 


ae 1 


r2 


at (=a) eee (1 —o) 


Equation (19) can be written in non-dimensional form as follows: 


dn, dn, @ [| (L+0)(1—20) 
a + ae [a exe 5) Lae 
ee te es woe = 0. (20) 


While equations (17) can be written under the following forms: 


a) Y, = 0 [displacement zero] or 7, =0, ) 
ee = 0 [total slope zero] or = 0, 


c) = 0 [slope due to bending zero] or 
dy, (1 + o) (1 —2 0) @ jdm 
des d—a) Mal eam se 
dX, (21) 


d) ah = 0 [moment zero] or 


d? 
wate on =O 


d 
e) - + X, = 0 [shear zero] or 


an, , B® (Le G) == 2:0) | eure 
dg + | f (1 — 0) a aa Us 


* In fact from equations (15) and (5) the following relationship is derived: 
L 


1B, L 
w? | [X\(x)]2 de + | [Y,(x)P af a Ret 5" | eae dx + c| | X,(x) + eal dx. 
0 0 : 


? See: E. Volterra, Ing.-Arch. 24 (1956) p. 317. 
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In equations (20) and (21) 


a <8 w? 9 L? 
ny te ? gE == Th > Pp = as > 
(22) 
z r2 r2 ( 
a=2(1+o)F5, OO =a: 


The quantities ao and 62 in equation (20) are proportional respectively to the shear flexibility and 
rotatory inertia of the vibrating bar. If the shear flexibility, rotatory inertia and the quantity 
o are neglected equation (20) reduces to the classical equation 


dy, ool? L 


dé exe 
By assuming the solution of equation (20) under the form 
Mm = Ce 
and substituting equations (23) into equation (20) the following characteristic equation is obtained: 
qs 4 ab (ql ss ia az a (1 PSS Fae) ares (24) 
Equation (24) can be written under the form 
(2 +2) (22) = 0 (25) 


the two roots A? and 3 being 


«  (l+o)(l1—20) 
ta d—o) 


/ yen ac) (9) 
i (a) 
2 meee - 


2 
e+ Spats Pert) (oe : 
o (i=) Ate (ia) _ &(1 +0) (1—20) 

2 1 2 ee (i==5) . 


From the first of equations (26) it appears that J? is a positive quantity. 
From equations (26) one obtains the following relationship: 


@ (1 + 0) (1—20) 


RA 


@? (Teg) (2 6) | 
Cae Se) $ (26) 


R=” 


he = (i=) (1 — B% a). (27) 
From equation (27) it descends that 
fie OO if Ag ee 1S (28) 
feet. 2x 1 (29 
The general solution of equation (20) will be 
n,(é) = C, sind, & + C, cos, € + C; sinh A, € + C, cosh Ay & (30) 


with the understanding that if conditions (29) instead of condition (28) is verified, the following 
substitution: - 

1, = +77, (31) 
should be made into equation (30). The frequency equations and modal functions will be now 
derived for the following cases: 

a) Simply Supported Bar. The boundary conditions are in this case: 
dX, 


at ereaOsae Yo S310.) io 8 
dX, 
ve Spear hee Grae 0) ae 0 


from which the frequency equation 

sin A, sinh/, = 0 
or ; Ai ee or Ase 0 (32) 
and modal functions 


Ha C, sink € (Bias 12,3) fe) (33) 
follow. 
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b) Cantilever Bar. The boundary conditions are in this case! 


at Or Yo UF GV 
dX. dY. 
at eile re sae oO 
from which the frequency equation 
Pe Dd? sin Ay, sinh Ag, _ 
2+ 2 ao a (« — on cos A,, cosh Ay, — ae (« + 6?) irae =) (34) 


(ete oer) 


and modal functions 


° . 2 
nan = Gu (Are a Maa (He + Fo) sin aes 


Aik Aok 
eee, a 1 %% 
—Ag, (it. — 9) sinh A,,;, : + ( +3 +34) cos A,, (35) 


+ ( 2 =e cosh Jay (cos A,,, € — cosh Ap, 5| 
follow. 


c) Built-in Bar. In this and in the following cases the length L will be regarded as half 
the length of the whole bar. The appropriate boundary conditions for symmetric or antisymmetric 
modes will be specified at the center (€=0) and at one end (£=1) of the bar. In the case of a 
bar built-in at both ends the following boundary conditions have to be verified for symmetric 
modes: 


ey. dY, = xe = 
at =0: =1=0, apie (36) 
at = 1; Y¥, =o, G0; 
From equation (35) the frequency equation 
2 2 
dns (A + #°) cos A,,, sinh A,, + Ag, (4. ae x) sin A,;, cosh A,, = 0 (37) 
(eae Sheen) 
and modal functions 
= GC, (cosh Ay, cos A,, & — cosd,, coshA,,&) (k =1, 2,3, ...) (38) 
follow. 
For antisymmetric modes the following boundary conditions have to be verified: 
at FO, a2 On eee 
dx (39) 
at =? G Ya =s0h, Xa 0r 
From equations (39) the frequency equation 
@? : ? , 
Aik (2 + ae “| sin A, ,, cosh Ay, & Aas (Be — #°) cos A, ,, sinh A,,, = 0 (40) 
and modal functions 
my, = C, (sinh A,,, sin A, ,, € — sin A,,, sinh A,, €) (41) 


follow. 


d) Free-Free Bar. For symmetric modes the following boundary conditions have to be 
verified: 


caye dY, 
t — © ou ee aes — 
a é=0: St=0, or ao 
42) 
dX dyY. ( 
t == , eres hs —- Sas = 
a é=1: St=0, BATE =x | 


} It is supposed in this case that the fixed end of the bar be situated at € = 0, while the free end at € = 1, 
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From equations (42) the frequency equation 


fan (le 5 be fo) 008 A ;, sinh A,, + Aan (Ae + «) sin Ay, cosh Le) (43) 
and modal functions 
M = G, (Ay, sinh Ag, cos Ay, & —Ag, sin A,;, cosh A,¢) (44) 
follow. 
For antisymmetric modes the following boundary conditions have to be verified: 
dX. 
Le = 0: = eee 
a € = Ss ide 0, 
4 (45) 
eo eke dY, n> | 
at &é=1: a Ue i 
From equations (45) the frequency equation 
Pd ses ae ee eo; 
(Ate — 34) sin A,, cosh A,;, —A,, (28. #° cos A,;, sinh /,, = 0 (45 a) 
and modal functions 
Ne = C, (Ay, c08 h Ag, sind, € + Ag, cos A,, sin h Ag, €) (46) 


follow. 


e) Bar Built-in at One End and Simply Supported at the Other. If L represents 
the whole length of the bar the frequency equation and modal functions for this case are the same 
as those of the antisymmetric modes of the bar with both ends built-in which have been dealt 
previously [see section c) of this paragraph]. The frequency equations and modal functions which 
have been derived apply to lateral vibrations in the three-dimensional case [equations (1), (2), 


(5), (6)]- 

In the case of lateral vibrations in two-dimensional elasticity [equations (10) and (12)] the 
same expressions for frequency equations and modal functions are valid, but the following changes 
have to be introduced: 


l—o 5 : 
( ) the constant has to be introduced in the 


1 
1. Instead of the constant Gao ==20) (1 — 0?) 


formulae. 


2. The roots A?, and /3, of the characteristic equation are given by the following formulae: 


a @? 
} 62 


Z = 
n= Gt —o9) + / otf —a—aaf + Fa—oy, 


1 ie: 2 


3. Lateral Forced Vibrations. By introducting into the variational principle the work done by 
the external forces the equations of motion are obtained in the case of forced vibrations. In the 
case of lateral vibrations equations (1) and (2) transform into the following equations: 


(l—o)E° , Py a » Pry 
(+ o)—20)" ox? per agatee ll o> (47) 

d?yy Oy =» vy ~——Py(x, t) 
Ce a = oe Tier 2 

(l—o)E 9 O'Ux = ee Ole 
C4 gG— 29) ta leo ae 2" ae (48) 
G PM 4G Abe vz P(X, t) 

axe! Bye @ Or A 


30 
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where p,(x,t) and p,(x,t) represent respectively the external forces per unit length of the bar 
acting in the direction of the y and z axis. A represent the cross-section area of the bar. Let 


Aes) = SX) ful)» wel) =F Xl) ful 


= es (49 
rylest) = 3 Youle) ful» v.C08) = 3 Vaal) fuel 7 
where the functions X,,(x), Y,,(x) are the solutions of the differential equations 
a aus at 7 6X, — 0! Tee aus 
69 4 621 pat ¥, 


which satisfy also the usual homogeneous boundary conditions to which the functions /,(x,t), 
vy (x, t) and w,(x, t) and v,(x, t) are subjected. Then, there exists an infinite sequence of charac- 
teristic values of w: 


Oy < Wy << Wg °° LO, <°** << W,°°° (0; > 0c as n>) 
to which correspond the non trivial solutions 
Xii(x) > Xyo(%), Xag(x)>-- +> Map (%). +--+, Xi a(x) ---> 
Yi (%)-+~— Yiolx) 2 Yag(*) <--> Yuy(x)> -+-> ¥q,@) =: 


Upon normalization of the eigenfunctions X,,(«), Y,,(«) which represent the normal modes of 
vibration of the system one has: 


ion 0 for hS2k", 
is [ [¥aul0) Yaw (2) + PX) Nelo] dx = (51) 
3 Wintorthwstkes 
Substituting equations (49) into equations (47) or (48) one reaches the equations 
< (—o)E d2 X41. (x dY. ~ df; p(t 
ae la +o) (120) Iie oa G X,,(x) — oe ih 4(t) = @ r De X4,(x) eee ’ 
a d?Y;4(x) dX 4(x) X fe a 
2 a a d? t mt 
ps [oe + 6 SH) 640 = o> BOE ay ose 
and in view of equations (50) equations (52) become 
ioe) ad t oe) 
> X@) SEO = — 2 Sot, Xi) ful» 
k= 1 1 (53) 
oe af; A(t at . 
Sy Yut) fu) _ Pe) Doh Yau) fix 


Multiplying the first of equations (53) by X,,(x)dx, integrating between 0 and L, multiplying 
the second by Y,,(x)dx and integrating between 0 and L adding subsequently the two equations 
together, by virtue of the orthogonality relations (51) one is lead to the following equations: 

L 


d2 
Fe + ot ful) = Ql) = sage | pl) Yaule) de. (54) 


0 
The problem of determining the response to the system of the external applied force is therefore 
reduced to that of finding the particular solution satisfying equation (54) and conforming to the 
initial conditions. In the same way the problem of longitudinal forced vibrations can be solved. 


fs Longitudinal Eigenvibrations. In the case of longitudinal vibrations expressed by equations (3) 
gaa as v0 1) = Kale) fC); 
A,(x,t) = Y,(x) f,(t) , (55) 
L(x, t) = Z,(x) fad) « 
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Substituting equations (55) into equations (3), one obtains 


d*Y,(x) , ox: dX,(x) 
Gr —(A-+ 26) Y, a 
mat 2) XG) A ALG) 1 at 
er? Y,(x) RO mde oo Re 
d? ao X,,(x) 
Gry —(A+26 ae oe? 
¥ ) Z,(x) —A—F A Y,(x) L af) _ , (56) 
+ Z,(x) iF WAGLE Sears BP 
Q+26) ae) im Re ee 
oe X,(x) ~ FO dee 


Since the constant «,is a real quantity, the longitudinal eigenvibrations are defined by the follow- 
ing equations which descend from equations (56) and the boundary conditions (7): 


SHO + of) = 0, 
pee Eee iap ars PHO) 4 Z(x) + who 1? Yx(x) = 0, 
G13 SHO) _ (4 4.2.6) Zy(x) —2 2) 2 Ye) + whe Ayla) = 0, sa 
(A +26) SX) 4 OO) 5 Fa) + apo Xx) =0, 
ae = =i 
ze] a=9 CEA 0x le (58) 
G pe hy za, 0 


The equations of motion (57) together with the boundary conditions (58) are satisfied for a triple 
set of discrete eigenfrequencies ~,, to each of which corresponds a fundamental mode given by 
the eigenfunctions X,,, Y,,, and Z,,. From please (57) and (58) the following orthogonality 
property for the eigenfunctions X,(x), Y,(x), Z.(«) can be derived :? 


L 
PX al?) anA®) tp 1s You(*) Yau(*) + 15 Zay(%) Zp,(x)] dx =0 for en F Man, (59) 
where X,,,(x), Y,(x), Zo,(x) are the eigenfunctions corresponding to the eigenvalue (,,. 


If the bar has a circular cross-section equations (3) and (7) reduce to the following equations 
of motion and boundary conditions?: 


cp th E “ oE dry yp Dy 
(ate (it oles) (1 25) a) dx 8 ae Bs 
(lL—o)E ary , 2vE dhe dry ) 


(+ 0)(1—20) de® * +0) —20) dx © da?’ 


1 In fact from equations (56) and (7) the following relationship is derived: 


L L L L 
0 03 ‘ | [Y,(x)]? dx + r} | [Z,(x)]? dx + | [X,(x)]? al = Cie f eae dx + Gri [ aah 
0 0 0 0 0 
x dX, (x) 
426 | ve) Z3(x) 4 ae ea | ¥,(x) + Zi) + 2 eo] a 
0 0 


2 See: E. Volterra, Ing.-Arch, 23 (1955) p. 410. 
3.0% 
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eee 
| EA 


a 7 me ess gesihc (61) 
E ae eae OP =O 

Let 
V (x,t) = X(x) f(t), A(x, t) = R(x) f(t) - (62) 


Substituting equations (62) into equations (60) and (61), one obtains 


oF + oF f) =0, 


dt® 
@2R(x) E cE dX(x) = 
Cerin = (ee Daiaea) in) ep asprey (aa) ne Cita aaa (63) 
(l—o)E  dX(x) 20E dR(x) » 
(cia ea) Sees) DN) 
and 
‘ze 0 
“ise 5 a ee O, ee lh. (64) 
oO 
E eas) R)X =0 


By eliminating R(x) between the last equations (63) the following equation in X(x) is obtained: 


d2X( wo? 
a 7” “aE X() ; } 


(p56) (Lana) CA aX) eae) Ea peels (65) 


\ Gr? d?xX Gr? dx 
(+9 0—20)(8)S A (lao 22) 0 


The first two terms of equation (65) are the classical terms, while the second and third lines represent 
the terms due to radial shear and to radial inertia respectively. 


Equation (65) can be written in non-dimensional form as follows: 


—fa eae cel (1 ss emma 


+7 ([O? — (1 +4) 1 — 20) OF04] = 0. 


(66) 


In equation (66) the terms containing 64/x are the terms due to radial shear, while the terms con- 
taining 6* are the terms due to radial inertia. If these terms are neglected equation (66) reduces 
to the classical form 


d2 
ee + On =0. 


Equations (64) expressing the boundary conditions can be written in non-dimensional form as 
follows: 


a)! X09) Or. 7 


Kd 


aa dy UP oy 20), 2 0? a|dyn _ 
b) R=0 o9 Ge lo (=e) =a) dire ec) wae 9 
hecho din pa (0) 20) 
ee ap MT RIM Grease (lo) eet (67) 
dX 20 2 dy (1 + o) (1 —2 0) a 1 a) dyn 
Oe ie (6) ae eae {| Tsao & Aree | tee? 
dX dy 


XXV. Band 1957 Volterra and Zachmanoglou: Free and Forced Vibrations of Straight Elastic Bars 433 


In equations (66) and (67) 


x rp IEP 
Geerage aurea! aoe 
Er? 2 r2 (68) 
mec gar Tas 


By assuming the solution of equation (66) under the form 


and substituting equation (69) into equation (66) the following characteristic equation is obtained: 


1 a a (1 + co) (1—2 0) 
MR Gaya Cola (1 —o) 
Pe _ gi ity G—20) «\* , ey 
|a—s) i (1 — 0) | = ’ 
Equation (70) can be written under the form 
(2? + Aj) (4? —23) = 0 (71) 
the two roots A? and j2 being: 
BR=—A+B, B=+A+B, (72) 
where: 
me! 1 4 ae je 0) 26) 
a 2 laa ot ® E C=.) I 


2 or oa os ol a (@+oe(—20) « oe 


ca 1 1 o 
: laa ue a! (16) 1—o) & do) 


From equations (72) it descends that 


: 1 (1 + o) (1— 20) pees 

f 2 : 2 
Teale (ly ee) a (73) 
then eS = 05 

“¢ rite. (et yl 20) pee ie oa ae) 2 A Dees 

Be eee G9 fe ca (74) 
then Aa OL eapees 0” 


The general solution of equation (66) will be 
n(é) = C, sind, € + C, cosd, & + C; sinh A, € + C, cosh, & (75) 


with the understanding that if condition (74) instead of condition (73) is verified, the following 
substitution: 
tag = 1A, (76) 
should be made into equation (75). 
The frequency equations and modal functions will be now derived for the following cases: 


a) Pinned end Bar. The boundary conditions are in this case 


dR 

at 2 se (Dp a — 10 a 
dR 

at aah Xxe—1()ie eo. 


from which the frequency equation 
sin/, sinh A, = 0 
ODneA = 6 tor. Ag 0 (77) 


and modal functions 


n, = C,, sink wx & (78) 
follow. 
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b) Cantilever Bar. 
The boundary conditions are in this case?: 


at E208 “EX 0 he 08 
dR dX om ae 
at Federer he os (Rao ae 


from which the frequency equation 


—[Q—o) Re + +o) G—20) 8) | (L—o) Mj, —(1+ 9) — 20) & 
oa a o) 12, = (1 + 0) (120) & | (io), + +c) (1—20) cos A,,, cosh Ap, 


- A us sind psinl 1450 a 
and modal functions 
Pears le ie an (A ae toh ! oi) sind, é 
Malte AGE atag fk AEG 
! ( ae kS a) cosh Ap,| (cos 44, — cosh Ay, »| } 


follow. 


c) Built-in Bar. . 

In this and in the following cases the length L will be regarded as half the lenth of the whole 
bar. The appropriate boundary conditions for symmetric or anti-symmetric modes will be specified 
at the center €=0 and at one end =1 of the bar. In the case of a bar built-in at both ends the 
following boundary conditions have to be verified for symmetric modes: 


dx dx 20 
at & z= ilo SS =O. R=0 


from which the frequency equation 


(1 + o) (1 — 20) 
d= oe* 


Ayi, Asi, + 


| cos A, ;, sinh A,, 
(81) 


1 1—2 
+ hop an ( Tes Gj g) vi sin J,, cosh A,, = 0 
and modal functions 
Fa Ne = CG), (cosh Ay, cos A, & — cos d,;, cosh A,, &) (82) 
ollow. 
For antisymmetric modes the following boundary conditions have to be verified: 
at | P= 0 0 ne 
dx 
at é == is ne — 0 > R == 0 
from which the frequency equation 
1 1—2 
Auk te + ( toe 5) 2) @j,| sin 1, ;, cosh Ay, 
1 1—2 (83) 
=a As ie — ( Tas >) g) D;, [eos Ark sinh Ask = 0 
and modal functions 
n, = C, (sinh dy, sind, § — sind, ,, sin h Ay, €) (84) 


follow. 


* It is supposed in this case that the fixed end of the bar be situated at € = 0, while the free end at € = 1, 
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d) Free-Free Bar. 
For symmetric modes the following boundary conditions have to be verified: 
ax i do 


at e= 0; ae ares Re ea ae 
dR dx 20 
at B= Ils me & ae (1—o) = ()) 


from which the frequency equation 


Jax 78 jl (eb o)( ie 


a ot sin A,,, cosh Ay, 


2k (l—o) 
(85) 
1] —2 
+A), a (l TSS °) p cos /,, sinh A,, = 0 

and modal functions 

my, = C, (Ay; sinh A,;, cos A, & —Aq;, sin A,,, cosh Ay, &) (86) 
follow. 

For antisymmetric modes the following boundary conditions have to be verified: 
at fe et 
dx 
dR dx 20 
at eee alles Sie ee (eee 
from which the frequency equation 
1+oj)1—206 ; 
Aik vA ( Hee x ) oi sin A,,, cosh Ay, 
(87) 
1+o)j)(1l—206 : 
Aan [ibe + § a ” i} cos Ay sinh Jay 

and modal functions 

Ny, = C,, (Ag), €08 Ay ;, sinh Ay, & + A,;, cosh Ag, sin A,),€) (88) 


follow. 


e) Bar built-in at one end and pinned at the other. 


If L represents the whole length of the bar the frequency equation and modal functions for 
this case are the same as those of the antisymmetric modes of the bar with both ends built-in 
which have been dealt with previously [see section c) of this paragraph]. The frequency equations 
and modal functions which have been derived apply to longitudinal vibrations of a bar of circular 
cross-section according to equations (60) and (61). 

In the case of longitudinal vibrations in two-dimensional elasticity [equations (11) and (13)] 
the same expressions for frequency equations and modal functions are valid, but the following 
changes have to be introduced: 


1. Instead of the constant 


a a 2 20) the constant i has to be introduced in the 


formulae. 


2. The roots 4? and 2 of the characteristic equation are given by the following formulae: 


ee ne ae By (89) 


where 


ib [ice 4 2 a a 
B=V7le—-2 |e t+a—o} | E p< (1—o)|. 
5. Torsional Vibrations. In the case of torsional vibrations expressed by equations (4) and (8), let 


by(x5t) = Tx) fA), — A,(% t) = (x) Ad) - (90) 
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Substituting equations (90) into equations (4), one obtains 


2 dt*(x) 
Cr: Jp G u(x) —& T(x) Sly POTS aie 
gr? Hx) OM a ” 
or} Te) fQ ae ~~ 


and since the constant jw, is a real quantity! the torsional eigenvibrations are defined by the follow- 
ing equations which are deduced from equations (91) and the corresponding boundary equations (8): 


d? f(t 
PHO + ob felt) =0, 
Gri 9) _ Giz) — CT) +0 lr x) =, (92) 
Ge AE 
G ry at T(x) —G t(x) + @ wiry T(x) =0 | 
dT 
gro 
ata 0 (93) 
0 
fs| 
The eigenfunctions T(x), t(x) satisfy the following orthogonality property: 
if 
f [re ta(x) t,(x) + 7} T(x) T,(x)] dx = 0 for. w2, 4 Wem. (94) 


In the particular case in which the cross-section of the bar is circular, the following conditions are 
satisfied?: 


Ty =T,=T, —fly=A,, Ux) = — T(x) 
and equations (92) and (93) reduce to the equations of the Classical Theory 


d? 
Te + off) = 0. 


aaa 
6 TO + ooh Te) =0, 


|= ate i — 1 bemeve—alors (96) 


(95) 


dx 


The material presented in the paper is based on an investigation conducted at the Rensselaer 
Polytechnic Institute, Troy, N. Y, under the sponsorship of the Office of Ordnance Research, 
Contract No, DA-30-115-ORD-709. Project No. 454.13. 


(Eingegangen am 16. Dezember 1956.) 


Anschriften der Verfasser: Prof. E. Volterra, Troy N. Y. (USA) 1567 Tibbits Ave; 
E. C. Zachmanoglou, Troy N.Y. (USA) 
1 In fact from equations (91) and (8) the following relationship is derived: 


L L L L L 
: 3 : 
nef | [T(x)}? det | [e(x)]? a} on ieee e+ 6H | ae dx + G | [T(x) + #(x)]? dx. 
0 6 ) i) 0 


2 See: E. Volterra, Ing.-Arch. 23 (1955) p. 410. 
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eit -winkel. _ Behandlung von ‘Rahmenwerken mit ein- und iachatachct Symmetrie. — Ein peesadatss Ver- 
fahren fiir eingeschossige Rilaeareda —_ Die Handhabung des Cros pVevishivems bei sréferen 


Stockwerkrahmen, Rahmentriger und bei Darehlauttrigern auf elastisch senkbaren Stiitzen. — Einflub von 
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Temperaturinderungen, Schwinden und Normalkriften. — EinfluBlinien. — Veriinderliche Trigheits- 


Daaientes — - Anbangs Tabeen IIx. Namen- oad Sachverzeichnis. 


‘ ’ 


: .. Die Verfasser bemithen sich, das von Cross entwickelte ataks Iterationsverfahren streng in seinem 
Bes Sinn darzustellen, auch fiir solche Fille, deren Lésung Cross nicht angab oder nur- andeutete. | Es kann 
_ festgestellt werden, daB es panes gelungen ist, in klarer Shetek, pseeeclidnd diesen Vorhaben ge- 
3 “a recht zu werden, und zwar sowohl fir Rahmentragwerke mit yesemeee unverschieblichen als auch fiir 
; ; | solche mit Vetsckiohliclien Knotenpunkten sowie fiir Durchlanfirager auf starren Stiitzen. Dazu wurden 
| | aueli Berechnungsverfahren angegeben, die das Crosssche Verfahren abkiirzen a dabei teilweise von 
4 ; ; Mitteln der Zusammenfassung mehrerer Iterationsschritte (Konvergenzbeschleunigung) Gebrauch machen, 
| s ance auf der Anwendung einfacher statischer Bvsibme beruhen, Mit den fiir Rahmen und Durchlanf- . 
trager auf starren Stiitzen gewonnenen Erkenntnissen werden auferdem Durchlauftrager auf elastisch 
ge Stiitzen und mit gegebenen Auflagersenkungen, der Einflu8 von Temperaturanderungen, des 
Schwindens und der Normalkrifte behandelt. Weiterbin erfahrt man, wie HinfluBlinien ermittelt und ver- 
anderliche Tragheitsmomente beriicksichtigt werden. 
"Die zablreichen Beispiele tragen zur Verstandlichkeit des Buches bei; sie zeigen aber auch die Grenzen des 


Cros s-Verfahrens, wodurch die Gefahr eines MiBbrauches verhindert wird. »V DI- Zeitschrift 
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Von Dr.-Ing. Wilhelm Fliigge, Professor of Engineering Mechanics, Stanford University. “Zweite, 
neubearbeitete Auflage. Mit 121 Abbildungen. VIII, 286 Seiten Gr.-8°. 1957. _ Ganzleinen DM D8 50c oe 
Inhaltsiibersicht: I. Allgemeine Grundlagen: Einfihrung. ‘Membranspannungszustand. Transforma- x 
tion der Langskrafte. Langskrafttransformation in schiefwinkligen Koordinaten. Transformation der 
Momente. Theorie der Beanspruchung von Stahlbetonschalen. — IL. Membrantheorie der Rotations- et 
schalen: Allgemeines. Drehsymmetrische Belastung. Fragen der Formgebung. Unsymmetrische Be- 
lastung. Forminderungen. Formanderungsarbeit. Statisch unbestimmte Systeme. — III. Membrantheorie a 
der Zylinderschalen: Rohre und Tonnendacher. Formanderungen. Statisch unbestimmte Zylinderschalen. 
Vieleckkuppeln. — IV. Membrantheorie fiir beliebige Schalenformen: Gleichgewichtsbedingungen fiir eine _ 
Schale beliebiger Form. Elliptisches Paraboloid. Hyperbolisches Paraboloid. Konoidschale. Der Spannungs- — 
zustand affiner Schalen. — V. Biegetheorie der Kreiszylinderschale: Die Differentialgleichungen der Zylinder- 
theorie. Lésung fiir die an den Randern x = const belastete Zylinderschale. Theorie der ‘Tonnendacher. aie 
Behiltertheorie. — VI. Biegetheorie der Drehschalen: Differentialgleichungen. Kugelschale konstante1 s 
Wandstiirke. Schale mit beliebiger Meridianform und veranderlicher Wandstarke. Kegelachale. — VII. Theo- 
rie der Faltwerke: Membrantheorie. Biegetheorie. — VIII. Die Stabilitat der Schalen: Grundlagen der Theo- 
rie der Schalenknickung. Stabilitat der Kreiszylinderschale. Stabilitat der Kugelschale, IX. Schwin- — | 
gungen drehsymmetrischer Schalen: Dehnungslose Schwingungen. Vollstandige Schwingungstheorie fiir die 
Kreiszylinderschale. — Literaturiibersicht. — Sachverzeichnis. ps oa 


- . ; . > i 
Das Buch behandelt die Membrantheorie der Drehschalen und der Zylinder, eine allgemeine Membrantheorie fir beliebige 
_ Schalenformen, die Biegetheorie der beiden wichtigsten Schalentypen, das Knicken oder Ausbeulen diinner Schalen und die 
Schalenschwingungen, Ein Abschnitt ttber die Spannungsberechnung in Stahlbetonschalen und ein AbriB der Faltwerktheorie 
werden dem Stahlbetonmann besonders willkommen sein. Die Darstellung benutzt die dem Zweck der Sache entsprechenden 
mathematischen Hilfsmittel, ist jedoch stets darauf bedacht, der Mathematik eine dienende Rolle zuzuweisen und sie nie. 
zum Selbstzweck werden zu lassen. Die zahlreichen Beispiele zeigen, daB sich alle dargebotene Theorie in praktisch anwendbare __ 
Ergebnisse umsetzen laBt. . j : 5 a ; yom os of 
Das Buch wendet sich an einen doppelten Leserkreis: an die am Weiterbau einer technischen Schalentheorie interessierten : 
Forscher und an die diese Theorie anwendenden Ingenieure aller Fachrichtungen. Theoretisch interessierten Studenten des In- : 
genieurbaus, des Maschinenbaus und des Flugzeugbaus wird das Buch ein Fihrer sein beim Eindringen in ein aktuelles Gebiet — 
der technischen Mechanik. __ ; / ee Sh a : 
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Elementare Schalenstatik = 
Von Dr.-Ing. habil. Alf Pfliiger, Professor an der Technischen Hochschule Hannover. Zwei te erweiterte | 
Auflage der ,,Einfiihrung in die Schalenstatik“*. Mit 56 Abbildungen. vu, 112 Seiten. Gr.-8°. 1957. 
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Inhaltsiibersicht: I. Einleitung. Allgemeines. Rechnungagrundlagen, — 1. Membrantheorie der Rota- E 

tionsschalen. Geometrie der Rotationsschalen im Kuppelbau. Rotationsschalen im Behalterbau. — : x 

Il. Biegetheorie der drehsymmetrisch belasteten Rotationsschalen. Gleichgewicht am Schalenelement. 

Elastizitatsgesetz fiir die SchnittgréBen. Kreiszylinderschale. Naherungslésung fiir beliebige Rotations. 

schalen. — IV. Membrantheorie der Zylinderschalen. Kreiszylinderschale mit -waagerechter Achse unter ee 

Eigengewicht. Zylinderschalen allgemeiner Form. — V. Membrantheorie allgemeiner_ Schalen. Bezeich- | fs z 

nungen und geometrische Bezichungen. Gleichgewichtsbedingungen. Differentialgleichungen. Hyperbo- 

lische Paraboloidschale unter Schneedruck. — VI. Einzelheiten des Spannungszustandes. Berechnung der 
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Spannungen aus den SchnittgréBen. Abhangigkeit des Spannungszustandes von der Schnittrichtung. — VII. oe 
Anhang. Schrifttumshinweise. Zusammenstellung von Lésungen der Membrantheorie. — Sachverzeichnis. — & 
‘ A . ! Loe ee 
Das Buch vermittelt fiir den Studenten und den in der Praxis tiitigen Ingenieur die Grundtatsachen der fiir den Statiker nicht _ “. of aa 
leichten Schalentheorie. Es liefert so eine Grundlage einerseits fiir das Verstindnis der schwierigen Spezialliteratur, anderer- a oe 
seits fiir die Lésung einfacher Aufgaben der Praxis. Die bewufte Beschrankung des Stoffes auf die ,,elementare“ Schalenstatik _ too as 
grenzt den Aufgabenkreis des Buches gegeniiber dem des oben angezeigten Werkes von W. Fliigge ab. Nach einerausfiihrlichen ate 
Erlauterung des Begriffes einer Schale und der Annahmen fir ihre Berechnung wird zuerst die Klasse der Rotationsschalen be- é th 
handelt. Aufer der einfachen Membrantheorie, deren Handhabung an einer Reihe wichtiger Beispiele des Kuppel- und Be- ea 
halterbaues gezeigt wird, gelangt auch die Biegetheorie wenigstens fiir drehsymmetrische Belastung zur Darstellung. Fir die ete 
Kreiszylinderschale wird diese Theorie in exakter Form, fiir die iibrigen Schalen in Gestalt der Geckelerschen Naherung ge- Di RE 
bracht. Es folgen die Membrantheorie beliebig gestalteter Zylinderschalen, die das praktisch besonders wichtige Gebiet der Se 
Tonnen umfassen, und die Pucherschen Gleichungen fiir Schalen ganz allgemeiner Form. Nach einigen Betrachtungen iiber die = a 
Bemessung von Schalen schlieSt das Buch mit einer Formelsammlung ab, in der fiir die Praxis eine Reihe immer wieder bendtigter Le ae 
Lésungen der Membrantheorie zusammengestellt ist. a ~ Se ia has me Se 9 vay 
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